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ABSTRACT 

The PdBI (Plateau de Bure Interferometer) Arcsecond Whirlpool Survey (PAWS) has mapped 
the molecular gas in the central ~ 9kpc of M51 in its ^^CO(l-O) line emission at cloud-scale 
resolution of ~40pc using both IRAM telescopes. We utilize this dataset to quantitatively 
characterize the relation of molecular gas (or CO emission) to other tracers of the interstellar 
medium (ISM), star formation and stellar populations of varying ages. Using 2-dimensional maps, 
a polar cross-correlation technique and pixel-by-pixel diagrams, we find: (a) that (as expected) 
the distribution of the molecular gas can be linked to different components of the gravitational 
potential, (b) evidence for a physical link between CO line emission and radio continuum that 
seems not to be caused by massive stars, but rather depend on the gas density, ( c) a close spatial 
relation between the PAH and molecular gas emission, but no predictive power of PAH emission 
for the molecular gas mass, (d) that the I-H color map is an excellent predictor of the distribution 
(and to a lesser degree the brightness) of CO emission, and (e) that the impact of massive (UV- 
intense) young star-forming regions on the bulk of the molecular gas in central ~9 kpc can not 
be significant due to a complex spatial relation between molecular gas and star-forming regions 
that ranges from co-spatial to spatially offset to absent. The last point, in particular, highlights 
the importance of galactic environment - and thus the underlying gravitational potential - for 
the distribution of molecular gas and star formation. 

Subject headings: galaxies: ISM — galaxies: individual NGC 5194) 
1. Introduction 
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Understanding the distribution and properties 
of the cold molecular gas is fundamental for de- 
veloping a coherent model for star formation as 
molecular gas is regarded the fuel out of which 
stars form. The need for a physical model on the 
few parsec to kilo-parsec scale has become more 
urgent as theoretical simulations and models can 
now incorporate a more complete handling of the 
gas properties on these scales (e.g. Robertson & 
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Kravtsov 2008; Gnedin et al. 2009; Christensen 
et al. 2012). In particular, Giant Molecular Clouds 
(GMCs) received tremendous attention as most 
massive star formation in the Milky Way is asso- 
ciated with them (e.g. Mooney & Solomon 1988; 
Mead et al. 1990). This clearly suggests a very 
close link between (massive) star formation and 
molecular gas which has also been demonstrated in 
recent studies of external galaxies to hold at kpc- 
resolution (Leroy et al. 2013; Bigiel et al. 2008, 
2011). However, at smaller scales it is not ex- 
actly clear what is causing or dominating the phys- 
ical link, small scale physics, like local turbulence 
within clouds, or large-scale process such as den- 
sity waves. 

Most of the knowledge in this area is based on 
our understanding of molecular line emission ob- 
served in the Milky Way where the role of galac- 
tic environment is harder to study than in nearby 
galaxies due to the distance ambiguity. In the past 
few years there has been mounting evidence that 
the properties of the molecular gas and GMCs in 
particular are not necessarily uniform in certain 
galactic environments (e.g. Bolatto et al. 2008; 
Gratier et al. 2012; Rosolowsky et al. 2003; Dono- 
van Meyer et al. 2011; Rebolledo et al. 2012; Hi- 
rota et al. 2011; Nieten et al. 2006; Espada et al. 
2012; Wei et al. 2012). However, none of these 
studies covered a large enough sample of GMCs or 
environments within a single galaxy to obtain ro- 
bust statistics and search for changes of the molec- 
ular gas or the environment that could quantify 
the deviations. 

As the most abundant molecule H2 cannot be 
directly observed, therefore the line emission of 
the CO molecule serves as the most common sub- 
stitute. Over the past decades several relations 
between CO line emission and line as well as con- 
tinuum emission at other wavelengths have been 
observed. While the underlying cause for these re- 
lations is often not well understood, studies that 
aim to test on which spatial scales such relations 
might exist, i.e. globally versus locally, bear the 
potential to unravel the physical processes under- 
lying the relations. Several relations between the 
CO emission and other tracers of the ISM as well 
as the stellar population have been established 
(mostly) on global scales, however, it is unclear if 
they persist on smaller spatial scales. Prominent 
examples are the relations between CO emission 



and (far-) infrared luminosity (e.g. Young & Scov- 
ille 1991), CO emission and radio continuum (e.g. 
Allen 1992), and last but not least the relation be- 
tween star formation surface density and neutral 
gas surface density, i.e. the Kennicutt-Schmidt 
relation (e.g. Kennicutt 1989, 1998; Bigiel et al. 
2008). The exact underlying causes for these rela- 
tions are not fully understood, and several physical 
processes could lead to such relations. Therefore, 
observations at high spatial resolution are criti- 
cal to shed light onto the underlying physical pro- 
cesses which might act at different physical scales 
and potentially allow for the disentanglement of 
cause and relation. 

In this paper we combine new observations of 
the CO (1-0) line emission at 1" resolution from 
the PdBI Arcsecond Whirlpool Survey (PAWS) 
with existing archival multi- wavelength data rang- 
ing from the radio to the UV to study the relation 
between molecular gas as seen in its CO emis- 
sion and other galactic components (gas, dust and 
stars). The paper is organized as follows. After an 
introduction of the PAWS survey (§2), we describe 
the datasets used and the steps taken to obtain a 
common astrometry in §3. In §4 we present our re- 
sults based on 2-dimensional maps, the polar cross 
correlation technique and a pixel-by-pixel analysis 
for the various ISM and stellar tracers used. The 
results are discussed in the context of processes 
acting and affecting the molecular gas on clouds 
scale (§5) and the implications for changes in the 
galactic environment are presented in §6. We sum- 
marize and conclude in §7. In appendix §A we 
provide a short summary of the relevant molecu- 
lar gas cooling and heating processes in relation 
to observations of the CO emission line. A de- 
scription of the polar cross correlation technique 
is given in Appendix §B. 

This paper is the first of a series describing re- 
sults from the PAWS survey. The PAWS data 
products will be publicly available from a ded- 
icated web-page^. Pety et al. (2013) find that 
about half of the CO emission arises from a diffuse 
component that is consistent with a thick molec- 
ular disk. In addition a detailed description of 
the observations as well as data reduction can be 
found there. Motivated by the results of this pa- 
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per that shows environmental changes in the rela- 
tions between the CO emission and other tracers, 
Hughes et al. (2013a) study the probability distri- 
bution functions of the molecular gas as a function 
of environment and find profound difference be- 
tween center, arm and inter- arm regions. A follow- 
up of this result is the work by Colombo et al. 
(2013) who identify about 1,500 individual GMCs 
in the PAWS data and find that several CMC 
properties show a strong dependence on galactic 
environment, and thus the large-scale dynamics. 
A comparative study of CMC properties in M51, 
M33 and the LMC by Hughes et al. (2013b) shows 
that the CMC scaling relations found are very sen- 
sitive to methodological biases and that distinct 
difference in the CMC properties between the 3 
galaxies studied exist. Finally, Meidt et al. (2013) 
provide an explanation for the varying CMC prop- 
erties found using dynamical pressure arguments 
and propose environmental dependences of their 
star forming capabilities. 

2. Motivation for the PdBI Arcsecond 
Whirlpool Survey (PAWS) 

While most stars in the universe have been 
born in galaxies with typical stellar masses of 
M^ = 10i0-6±0-4 Mo (e.g. Karim et al. 2011), 
detailed studies of the molecular ISM and GMCs 
- the entities out of which stars form - have fo- 
cussed on lower stellar mass systems in the Local 
Group. Studies of more massive spiral galaxies of- 
ten lack resolution and/or areal coverage to iden- 
tify statistically significant numbers of GMCs in 
different galactic environments (e.g. Bolatto et al. 
2008; Gratier et al. 2012; Rosolowsky et al. 2003; 
Donovan Meyer et al. 2011; Rebolledo et al. 2012; 
Hirota et al. 2011; Nieten et al. 2006). Thus high 
spatial and spectral imaging of (a significant part 
of) the molecular gas reservoir of a typical star 
forming galaxy has not yet been undertaken. 

We chose M51 ("The Whirlpool Galaxy") as 
the target of our survey as it is one of the clos- 
est {D - 7.6 Mpc, Ciardullo et al. 2002), face-on 
{i - 22°, Tully 1974, Colombo et al. in prep.), 
grand-design spirals and one of the most-observed 
galaxies in the sky resulting in an excellent cov- 
erage across the full electromagnetic spectrum. 
A summary of the basic parameters assumed for 
M51 is given in Tab. 1. Over 1,000 individual 



star clusters have been observed and character- 
ized based on HST imaging (e.g. Bastian et al. 
2005; Scheepmaker et al. 2007, 2009; Haas et al. 
2008; Hwang & Lee 2008; Chandar et al. 2011) 
and deep, high-resolution imaging of the galaxy 
in the radio continuum (Dumas et al. 2011), mid- 
to far-IR (Calzetti et al. 2005; Mentuch Cooper 
et al. 2012), Pa/Ha (Scoville et al. 2001), and 
UV bands (Calzetti et al. 2005) trace young stel- 
lar populations with a wide range of age sensitivi- 
ties (e.g., Calzetti et al. 2005; Murphy et al. 2011). 
In contrast to Local Group galaxies with resolved 
CMC populations, the molecular gas clearly dom- 
inates the ISM over the 9 kpc region PAWS covers 
(e.g., Schuster et al. 2007; Leroy et al. 2008) and 
the molecular gas is concentrated into spectacular 
spiral arms (e.g., Garcia-Burillo et al. 1993; Heifer 
et al. 2003; Schuster et al. 2007). This offers the 
unique opportunity to relate the physical proper- 
ties of the molecular gas and its clouds to differ- 
ent galactic environments, especially regarding the 
galaxy dynamics and radiation field. 

The fundamental goal of this large IRAM pro- 
gram is to resolve the CO emission from M51 
on the scales of individual GMCs. This enables 
a new approach to understanding star formation 
in galaxies, one centered on studying the phys- 
ical structures of interest rather than averaging 
over large areas. The inferred total molecular 
gas mass in NGC5194, the large spiral galaxy, is 
Mgas = 6.2 X 10^ M0. About 60% of this gas 
mass - comparable to the total molecular mass of 
the Milky Way - resides in the PAWS field-of-view 
(FoV) (Pety et al. 2013). Together with the exist- 
ing mult i- wavelength data, PAWS was designed 
to provide the data to study the evolutionary se- 
quence for GMCs, place strong constraints on the 
processes that drive cloud and star formation in 
disk galaxies, reveal the physical underpinnings of 
one of the most commonly used scaling relations in 
extragalactic astronomy (the Schmidt-Kennicutt 
'law'; Kennicutt 1989, 1998), and add a crucial 
data point regarding the lifetime of giant molecu- 
lar clouds, a lingering unknown that is critical to 
understanding the gas-star cycle in galaxies. 

These goals resulted in the following considera- 
tions for the observations with the IRAM Plateau 
de Bure Interferometer (PdBI) and the 30m single 
dish telescope: 

Sensitivity: The goal of identifying individual 



GMCs in M 51 determined the sensitivity required. 
In the Milky Way, most H2 is found in clouds with 

masses in the range 10^ 10^ Mq (e.g. Sanders 

et al. 1985; Blitz & Wilhams 1999). In order to 
detect GMCs with 1 x 10^ Mq at 5cr, a la point 
mass sensitivity ^ 2 x 10^ Mq is sufficient to re- 
cover clouds across this entire range. 

Resolution: The typical size of a Milky Way 
GMC is ^ 40 pc (Solomon et al. 1987) and thus a 
comparable spatial resolution is required to iden- 
tify and distinguish individual clouds. This reso- 
lution also marginally resolves the dust lanes ob- 
served in the optical in M51 (width ~ 1'^). 

Area surveyed: The survey area was chosen to 
incorporate three distinct environmental regions, 
the center encompassing the bulge plus the star- 
bursting ring, the spiral arms and the (relatively 
quiescent) inter-arm region. A substantial area 
in each regime is necessary to build up sufficient 
statistics on the molecular gas and its GMC pop- 
ulations. 

These considerations resulted in imaging of the 
^^CO(l-O) line emission from the central ~9kpc 
of M51 between August 2009 and June 2010 us- 
ing both IRAM facilities. The global properties 
adopted for M51 throughout the PAWS survey 
are summarized in Tab. 1. 

3. Data 

In addition to our own new CO (1-0) observa- 
tions, we made use of archival data that trace 
different phases of the interstellar medium and 
the stellar populations as well as star formation 
present in the disk of M 51a. Below we describe in 
detail where the datasets used come from and how 
they were aligned to a single absolute astrometric 
frame. A summary of the mult i- A data used here 
is provided in Tab. 2. 

3.1. IRAM PdBI+30m PAWS data 

We obtained a map of the molecular gas distri- 
bution in the central 11 x 7kpc of M51 as part 
of PAWS (PdBI Arcsecond Whirlpool Survey, PI 
Schinnerer). This IRAM Large Program observed 
a 60-pointing mosaic with the Plateau de Bure 
interferometer (PdBI) in all configurations and 
mapped the full galaxy system with the 30m single 
dish telescope in the ^^C0(1 — 0) line. The result- 
ing data cube of the combined PdBI and 30m data 



(i.e. short spacing corrected) has a resolution of 
1.16'' X 0.97'' (PA 73°) with a mean rms of 0.4 K 
per 5 kms~^ wide channel. A detailed description 
of the data acquisition and reduction is presented 
by Pety et al. (2013). We also use the moment 
maps derived from this PAWS datacube (for de- 
tails see Colombo et al. in prep.). The conversion 
from integrated brightness temperature per beam 
into molecular hydrogen (H2) gas surface density 
per 0.3" pixel assumed a Galactic conversion fac- 
tor of Xco = 2 X 1020(.^-2j^-ii^^-ig 

The PdBI observations consisted of two 30- 
pointing mosaics that were independently ob- 
served between August 2009 and April 2010. Each 
pointing within a mosaic was observed between 
two calibrator cycles ensuring good and similar 
uv coverage for each pointing as well as the two 
mosaics. Standard observations of quasars were 
used for the passband and atmospheric calibra- 
tion which was done using the GILDAS software 
package (Pety 2005). The 3mm receivers were 
tuned to 115.090 GHz corresponding to the red- 
shifted ^^C0(1 — 0) emission line with a LSR 
velocity of 471.3 kms~^. The 30m single dish 
was used in May 2010 to map a ~60 square ar- 
cminute field that contains the entire M51 system. 
Seven separate regions were observed in position- 
switch on-the-fiy mode; four of those cover the 
central 300" x 300", while the remaining three 
extend this map to the companion and the outer 
molecular spiral arms. The data was reduced in a 
standard manner and re-gridded onto a data cube 
with a pixel size of 4". The final PAWS data cube 
is obtained from a joint deconvolution of the PdBI 
and 30m datasets. It has 120 channels of 5 kms~^ 
width and a pixel size of 0.3 "/px. 

For the analysis presented here we use a 0th mo- 
ment map that was derived using a combination 
of the dilated mask technique and HI priors (e.g. 
Appendix B of Pety et al. 2013). A detailed de- 
scription of the moment calculation is presented by 
Colombo et al. (in prep.). Briefiy, a dilated mask 
is created by identifying emission above a given 
threshold of ta. Starting from these peaks the 
identified islands of emission are expanded to in- 
clude fainter emission above pa. If the central ve- 
locity of the islands was more than 30 kms~^ away 
from the central HI velocity (based on the 6" reso- 
lution robust HI cube from THINGS, see §3.6) the 
islands were excluded from the final mask. Due to 



the high angular resolution the PAWS CO(l-O) V' 
data samples much better the rising part of the 
rotation curve than the HI data causing a lot of 
real emission to be dropped. Therefore two masks 
were necessary to capture significant emission in 
the disk and the center, i.e. we used (t,p) of (4,1) 
and (10,1.5) for the disk and center, respectively, 
and had an HI prior for the disk and none for the 
center. 

3.2. HST archival data 

To compare our view of the molecular gas with 
M51's appearance in the optical and near-infrared, 
we use HST imaging that covers the PAWS FoV. 
The HST ACS Heritage (GO-10452; PI Beckwith) 
final reduced and drizzled mosaics (Mutchler et al. 
2005) in the F435W (B band), F555W (V band), 
F658N (Ha narrow band), and F814W (I band) fil- 
ters were obtained from the dedicated web-page^. 
No further processing was applied to the combined 
mosaics. To obtain a stellar continuum corrected 
map of the Ha (plus Nil) line emission we fol- 
lowed the approach of Gutierrez et al. (2011) by 
subtracting a linear combination of F555W and 
F814W scaled by 0.0878 from the F658N narrow 
line filter that contains the Ha line. Note that 
this leads to an over-subtraction in the central 30^' 
which is likely caused by the stellar bulge having 
a different spectral energy distribution than the 
disk. 

An HST I — H color map was obtained from the 
ACS Heritage observations and dedicated NIC- 
MOS H band observations in the F160W filter 
(GO-10501; PI Chandar). The NICMOS obser- 
vations were a set of observations at 18 positions. 
We dithered between two locations at each posi- 
tion. The NICMOS data was reduced in a stan- 
dard manner. Because the NICMOS observations 
were made using three different guide stars, their 
relative astrometry is not consistent. We first cor- 
rected for the rotation angle offsets between visits 
using stars that were common between the ACS 
image and the NICMOS images. We then solved 
for the positional offset between the NICMOS im- 
ages and the ACS mosaic by cross-correlating un- 
sharp masked versions of the NICMOS images 
with the ACS images. We applied these deter- 
mined offsets to each NICMOS image and used the 



standard DRIZZLE process to combine the images 
to form a mosaic. We formed an F814W-F160W 
(/ — H) color map by first smoothing the ACS 
F814W image to the resolution of the NICMOS 
F160W image and then dividing the F814W im- 
age by the F160W image. 

3.3. Spitzer archival data 

M51 was imaged from the near- to far-infrared 
by all Spitzer instruments as part of the SINGS 
Legacy project (Kennicutt et al. 2003). The 
Spitzer IRAC 3.6 /im and 4.5 /im images were re- 
processed by the S^G project (Sheth et al. 2010). 
Meidt et al. (2012) applied and developed a sepa- 
ration method of the stellar and non-stellar com- 
ponents present in these IRAC images based on 
the Independent Component Analysis that re- 
quires two adjacent IRAC bands as the only input. 
We utilized both the map of non-stellar emission 
at 3.6 /im and the contaminant- free old stellar light 
from Meidt et al. (2012). We used these maps of 
the old stellar light to subtract the stellar contri- 
bution to the 8/im image from the SINGS DR5^. 

To trace the recent star formation we use the 
Spitzer MIPS 24/im image from SINGS DR5 pro- 
cessed with the HiRes algorithm as described in 
detail by Dumas et al. (2011) resulting in an in- 
creased angular resolution of ~ 2A". 

3.4. Herschel archival data 

M 51 has been imaged by PACS onboard of Her- 
schel as part of a GTO program (PI Wilson, see 
also Mentuch Cooper et al. 2012). The archival 
photometric PACS data at 70 and 160 /im were 
reduced to level one using HIPE v6.0, and then 
used SCANAMORPHOS vl2.0 (Roussel 2012) to 
produce the final maps. As HIPE v6.0 was used, 
we converted the PACS images from fiight model 
(FM), 5 to FM, 6 by dividing by the factors listed 
in the PACS Photometer Point Source Flux Cali- 
bration Report vl.O, which are 1.119 (70/im) and 
1.174 (160/im). All images were converted to 
MJysr"^. The mean full width at half-maximum 
(FWHM) of the PACS point response functions 
are - 5.6^^(70 /im), and - 11.4^^(160 /im). The ac- 
tual images use pixels sizes of l'^ 

The [CII] emission line has been imaged as well 
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as part of this Herschel GTO program (PI Wil- 
son; Parkin et al. subm.). The version used 
here was reduced using the KINGFISH pipehne 
for PACS spectroscopy (Kennicutt et al. 2011) in 
HIPE V9.0.3063 using the calibration files in PACS 
calVersion 41. This is similar to the standard un- 
chopped pipeline, but includes some modifications 
for improved handling of transients in the data. 
In addition to the calibration steps described in 
detail by Croxall et al. (2012) the latest pointing 
calibration and the newest drizzle algorithms were 
used to project these data using HIPE 10.0.2538. 
The final angular resolution is ~11.5'^ 

3.5. GALEX archival data 

The FUV and NUV images are from GI pro- 
gram GI3-050 (Bigiel et al. 2010, and kindly pro- 
vided to us by PI Bigiel). M51 was observed for 
~10ks and details of the data reduction are pro- 
vided in Bigiel et al. (2010). GALEX simulta- 
neously images in the far-UV (1350-1 750 A) and 
near-UV (1750-2800 A) with an angular resolution 
(FWHM) of 4.0'' and 5.6'', respectively, and a field 
of view of 1.25° (Morrissey et al. 2005). 

3.6. THINGS HI data 

For comparison to the atomic gas we use the ro- 
bust weighted intensity map from THINGS (The 
HI Nearby Galaxy Survey; Walter et al. 2008). 
The intensity map has an angular resolution of 
5.78" X 5.56" (PA -68.0) and fully covers M51a. 

3.7. Absolute Astrometry and Alignment 

Unlike radio-interferometric images that natu- 
rally come with an absolute astrometry, UV/optical/IR 
images need to be tied to an external astrometric 
frame. Therefore we checked and, if necessary, 
corrected their astrometry as detailed below. A 
summary of our corrections applied is provided in 
Tab. 2. 

We used the VLA 20 cm image of Dumas et al. 
(2011) as our absolute astrometric reference since 
it encompasses the entire M 51 galaxy system with 
its angular extent of ^ 10'. Comparison to the 
2MASS J band image and a montage of the SDSS 
i band tiles (6 cover the full M51 system and are 
combined using the IRS A Montage tool^) obtained 



from the respective archives, suggests that the 
agreement of (extragalactic) point sources within 
the M51 system and around it is good. We 
find small non-systematic offsets between the J 
band image and the i band montage and there- 
fore decided to use the J band image as our as- 
trometric reference for the optical/NIR images 
since it also has the same geometric projection as 
the radio data. The final absolute astrometry in 
2MASS is < 0.1 " depending on the source bright- 
ness^. Before deriving the absolute astrometry we 
changed the projection of all images to the stan- 
dard radio interferometric system, i.e. CD matrix 
with rotation into CRPIX/CDELT/CRVAL de- 
scription without rotation, using the GIPSY task 
f itsreproj (van der Hulst et al. 1992; Vogelaar 
& Terlouw 2001). 

As some of the stars present in the images are 
moving during the period over which the different 
datasets were obtained, we used a combination of 
point sources (Galactic stars and stars/unresolved 
clusters within M51) using WCSTools (Mink 
1999)^ and the overall geometry of NGC 5194 and 
NGC 5195 to obtain our best alignment across the 
different wavelengths. We estimate that our final 
astrometry is good to about 0.15" for the high 
resolution HST images and slightly larger for the 
lower resolution IR images. In any case, we do 
not expect any profound systematic effects in our 
analysis (particular those done at 3.0" resolution) 
due to this small error. A detailed discussion of 
the astrometric accuracy of the PAWS CO (1-0) 
data is given by Pety et al. (2013). 

We verified that the relative astrometry of the 
HST images was correct. Comparison to the 
2MASS J band image showed that the HST Her- 
itage I band image (convolved to the same reso- 
lution as the J band image) is offset from the J 
band image in both RA and Dec. direction by 
0.1" and -0.4", respectively. We verified this de- 
termined offset by comparing the convolved HST 
ACS I band image to the SDSS i band montage. 
The same offset was applied to all other HST ACS 
images as well as the HST NICMOS H band image 
that had its relative astrometry tied to the ACS I 
band image. Note that Chandar et al. (2011) find 
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a similar offset from a comparison to Chandra X- 
ray imaging. 

The absolute astrometry of the IRAC images 
was found to agree very well with the 2MASS J 
band image, thus no correction was applied. We 
used a convolved version of the IRAC 8 fim. image 
to test the absolute astrometry of the MIPS 24 /im 
image, and find a small offset of 0.25^' and 0.5'' in 
RA and Decl. that we corrected. The same off- 
set was applied to the HiRes version of the MIPS 
24 /im image. Iteratively, we checked the absolute 
astrometry of MIPS 70 /im and MIPS 160 /im im- 
ages using the corrected MIPS 24 /im and MIPS 
70 /im images as reference, respectively. Again, a 
small offset of 0.375'' in Declination was found and 
corrected. 

For the Herschel PACS 70 /im and 160 /im we 
determined the offset based on the corrected MIPS 
24 /im image (plus MIPS 70 /im image for the 
PACS 160 /im map) and found small offsets of - 
2.5" and -2.0" that we corrected. The absolute 
astrometric uncertainty for the [CII] line map is 
larger (^ 1") as the centering was achieved us- 
ing both the MIPS 24 /im image as well as the 
SINGS Ha image convolved to 3". However, obvi- 
ous common structures can not be easily identified 
between these maps leading to the higher uncer- 
tainty. 

The astrometry of the GALEX NUV and FUV 
image was checked against a corrected and con- 
volved version of the IRAC 8/im as the overall 
morphology of these tracers is very similar. In 
particular, coincident point sources were utilized 
to determine an offset of 0.225" and 0.525" in RA 
and Declination that was corrected. 

4. Comparison of CO(l-O) to ISM and 
Stellar Tracers 

The study of potential relationships between 
molecular gas line emission, namely CO, and other 
tracers of the ISM and stellar populations from 
various ages at cloud scales can provide new in- 
sights in the underlying physical processes link- 
ing these emissions. In particular, the wide range 
of galactic environments (i.e. bulge/disk, spiral 
arm/inter-arm) sampled by the PAWS FoV could 
help to discriminate between different origins for 
existing correlations. After a description of the 
CO emission in the central ~9kpc as seen by 



PAWS (§4.1), we evaluate the correspondence be- 
tween CO emission and ISM (§4.2) and stellar 
(§4.3) tracers in this section. In the following sec- 
tion §5 we discuss potential physical causes for the 
relations and deviations seen. All the findings are 
placed in the wider context of galactic environ- 
ment in section §6. 

In the following we compare the distribution 
and intensity of the CO (1-0) molecular line emis- 
sion to other prominent ISM tracers such as 
the atomic gas emission lines, optical extinction, 
hot dust and polycyclic aromatic hydrocarbons 
(PAHs) emission, and radio continuum. In order 
to study the relation of the molecular gas to recent 
and former star formation we compare the distri- 
bution of the CO (1-0) line emission to line emis- 
sion and thermal radio continuum of HII regions, 
the light from (young) stellar clusters, the ultravi- 
olet continuum from young to intermediate aged 
stars (~ 100 Myr), as well as light emerging from 
the old stellar population. In order to capture 
the different aspects of the relationship between 
these tracers we present the data in classical 2- 
dimensional maps using the images at their native 
resolutions, a polar projection including a cross- 
correlation analysis, and a pixel- by-pixel analysis 
after convolving all data to a common resolution 
of 3.0". 

The polar coordinate system is very powerful 
for studies of spiral arms as it allows for an easy 
assessment of their properties. We use a linear 
47r presentation of the azimuthal angle (with '0' 
corresponding to west and a counter-clockwise di- 
rection) and a linear projection for the radius. Be- 
fore re-projection all images were convolved to the 
resolution of the 3.0" PAWS data cube and re- 
gridded to its pixel scale. We converted all fiux 
scales to units of "per pixel area" . The polar cross- 
correlation allows for a quantitative assessment of 
azimuthal offsets between two different tracers (see 
Appendix B for details). With this definition a 
$ > offset for, e.g., star formation tracers implies 
that gas is fiowing through the spiral mode inside 
co-rotation for a given perturbation to the gravi- 
tational potential. Several patterns (or perturba- 
tions) have been identified in the disk of M51 with 
different pattern speeds (e.g. Meidt et al. 2008; 
Vogel et al. 1993; Elmegreen et al. 1989; Tully 
1974). A detailed analysis of the various patterns 
and their associated pattern speeds (see appendix 



of Meidt et al. 2013) suggests that at any given 
radius (probed by the PAWS FoV) one is inside 
a corotation resonance and thus the expected re- 
versal of gas flow beyond the corotation resonance 
expected for a single pattern is not observed. 

To illustrate the differences between the differ- 
ent tracers and the CO emission more quantita- 
tively, we utilize pixel-by-pixel diagrams plus an 
additional environmental mask separating the cen- 
tral disk from the spiral arms and the inter-arm 
region (Fig. 1, bottom left). The density of the 
data points is presented on a logarithmic scale. 
Where appropriate we fit a linear relation using 
the IDL procedure ROBUST_LINEFIT, that per- 
forms a best straight-line fit accounting the errors 
in both coordinates as weights for the data-point. 
Lacking error maps for the different tracers we use 
the same weight of unity (rather than an error) for 
all data points/pixels. 

4.1. The Cloud-Scale Molecular Gas Dis- 
tribution in M 51 as seen via CO (1-0) 

The CO (1-0) line distribution exhibits two 
prominent spiral arms with the northern arm be- 
ing on average fainter than the southern arm (Fig. 
l^top). After winding for about 180 ° both gas spi- 
ral arms become significantly less prominent and 
even indistinguishable from the inter-arm in the 
case of the northern arm. In addition, the south- 
ern arm is showing a bifurcation exactly south of 
the nucleus that persists for larger radii. On the 
other hand the northern spiral arm shows a pro- 
nounced 'dent' north-east of the nucleus. Taken 
together this strongly suggests that the spiral arms 
do not represent a single pattern. The two spiral 
arms emerge from a ring-like distribution close to 
the center. The disk inside this ring shows fainter 
emission distributed in abundant smaller struc- 
tures. CO emission is present between the spiral 
arms, however, at an even fainter level and in a 
more flocculent fashion. 

In the polar representation, the CO (1-0) spi- 
ral arms are almost straight lines as expected for 
a log-normal spiral (Fig. 4, top left). However, 
close inspections reveals three breaks with differ- 
ing prominence for both arms at radii of ~27", 
~55" and ~90" with observed pitch angles of 
about 40°, 30°, and 10°, respectively. In addition, 
an azimuthal offset along the spiral arm at r ?^ 35^' 
is most evident in the northern spiral arm. The 



overall slope of the two spiral arms is not identical 
suggesting that the underlying gravitational po- 
tential can not be described by a pure m=2 mode. 
The molecular ring appears as an almost horizon- 
tal line for radii of 15^' < r < 20''. Below a radius 
of ^10" the sampling in azimuth is becoming too 
sparse in our presentation of the data for any firm 
conclusions. The difference in spiral arm pitch an- 
gle has already been noticed before by Henry et al. 
(2003) who analyzed the BIMA SoNG data (Heifer 
et al. 2003) and reported 'break radii' of 27" and 
60". The latter break at 60" is close to our in- 
ferred break at 55" . The discontinuity at 35" , on 
the other hand, appears strongly only at the high 
resolution of our data and is absent at the low res- 
olution of the BIMA SoNG map (e.g. Fig. 1 and 
2 of Henry et al. 2003). It is interesting to note 
that the radial profile of the azimuthally averaged 
torques (see Fig. 1 of Meidt et al. 2013) changes 
direction at radii of 20" (corotation resonance of 
the nuclear bar), 35" (start of spiral pattern), 55" 
(corotation resonance of spiral pattern) and 90" 
(start of second spiral pattern) consistent with the 
breaks observed in the CO distribution at 35", 55" 
and 90". Using a banana wavelet analysis Patri- 
keev et al. (2006) find values similar to ours for 
the pitch angles of the CO spiral arms and show 
that the pitch angles of the two spiral arms dif- 
fer by on average 10 ° over the PAWS FoV. Henry 
et al. (2003) suggested the presence of an addi- 
tional m=3 mode in the central 6kpc that could 
explain both the angular offset between the two 
gaseous arms and their asymmetry in intensity. 

4.2. Comparison of CO(l-O) to other ISM 
Tracers 

Two other tracers of the neutral gas can be 
compared to the distribution of the molecular gas 
at angular resolutions of ~ 6" , the atomic gas in 
its HI line emission and from the strongest cooling 
line of the neutral gas, [CII] at 158 fim. as imaged 
by PACS. Three different types of tracers of the 
dust in a galaxy have sufficient angular resolution 
for a detailed comparison to the CO maps: (a) 
an optical extinction map from HST, (b) maps 
of PAH emission from Spitzer, and (c) the hot 
dust heated by young stars. We expect varying 
degrees of correlation between the molecular gas 
distribution and these dust maps, given that exci- 
tation conditions likely vary across our FoV. While 



the cold dust is best traced at long infrared wave- 
lengths (e.g. Dale et al. 2012; Aniano et al. 2012), 
available maps have a fairly low angular resolution 
of, e.g., >12" for 210 /im onward for the Herschel 
SPIRE instrument. The non-thermal radio contin- 
uum is probing very different components related 
to the ISM and it can have several origins (SN, 
stellar remnants, radio jet, magnetic field etc.). 

4.2.1. Neutral Gas 

The strongest far-IR cooling line of the cold 
ISM in nearby galaxies is [CII] at 158 /im which 
can contain up to 1% of the total far-IR lumi- 
nosity (Stacey et al. 1991; Brauher et al. 2008; 
Gracia-Carpio et al. 2011). Typically most of the 
[CII] emission in nearby galaxies arises from the 
surface layer of molecular clouds where the gas is 
warm, dense and photo-dissociated by FUV emis- 
sion from either nearby hot (OB) stars and/or the 
interstellar radiation field (see, e.g., introduction 
of Stacey et al. 2010). Comparison to the CO(l-O) 
emission (Fig. 2, top left) shows a close correspon- 
dence between the peaks of both emission lines in 
the molecular ring out to a radius of ~ 20'^ except 
for a region in the north-east segment of the [CII] 
ring that has no counterpart in the CO (1-0) emis- 
sion. Between the ring and the spiral arms the 
[CII] emission exhibits gaps, similar to the ones 
present in the host dust emission. Along the spi- 
ral arms the most intense [CII] emission is patchy 
and is offset to the convex side of the molecular 
gas spiral arms. 

The robustly weighted intensity map of the 
THINGS data that images the distribution of the 
atomic gas shows very little pronounced structure 
in the central region of M51; a few peaks coin- 
cide with the molecular gas ring (Fig. 1, bottom 
right) . Along the molecular spiral arms the HI line 
emission becomes stronger reaching its maximum 
(within the PAWS FoV) roughly at the same loca- 
tions where hot dust emission is observed. Taken 
together the distributions of [CII] and HI strongly 
suggest that along the spiral arms both emission 
lines trace the dissociation process of molecular 
clouds caused by the massive young star forma- 
tion. Within the molecular ring the situation is 
more complex, as the [CII] emission is strongest, 
but no substantial HI emission is seen there. 



4.2.2. Optical extinction 

According to results for local Galactic clouds 
(e.g. Dickman 1978; Lada et al. 2007), maps of 
(optical) extinction are a strong proxy for the cold 
dust distribution. Comparison to the HST I — H 
image (Fig. 2, top right) shows that the optical 
extinction is a very good predictor of the location 
of strong CO (1-0) emission and its underlying de- 
tailed fine-structured morphology. This is most 
obvious in the spiral arms. In the inter- arm re- 
gion the correlation is less pronounced which in 
part is due to the presence of stellar clusters along 
spurs (see La Vigne et al. 2006; Chandar et al. 
2011) that are especially prominent in the bluer 
bands observed with HST (see §4.3.3). 

The polar cross-correlation of these two trac- 
ers shows no azimuthal offset nor any radial trend 
(Fig. 5, top right), consistent with both tracers 
being spatially coincident. The prominent (dust) 
clump in the inter- arm at an azimuthal angle of 
^240° and a radius of ~39" has a CO counter- 
part at a lower level. Inspection of the / and H 
band images suggest that the color of this clump 
is affected by a group of star clusters causing this 
apparent inconsistency. 

Overall, the molecular gas distribution above 
a H2 surface density Smh ^^^^^ of ~ 90 — 
100 M0 pc~^ is well predicted by the I — H color 
map. Detailed comparison on a pixel-by-pixel 
basis at 3.0^' resolution shows that the relation 
between the HST I — H color and CO integrated 
intensity is linear before it appears to saturate 
(Fig. 6, top row). While this trend is present in 
all environments, the slope in the central region is 
significantly shallower suggesting that the param- 
eters controlling the relation are changing. The 
similar slope for the arm and inter- arm regions 
means that the gas-to-dust ratio is not dramati- 
cally changing between arm and inter-arm. 

4-2.3. (Hot) dust emission 

The distribution of the emission from hot dust 
probed at 24 /im and 70 /im by the Spitzer MIPS 
and Herschel PACS instruments shows a bright 
central disk encompassing the molecular gas ring 
(Fig. 2, middle row). In addition, several bright 
peaks are detected along the molecular spiral arms 
with a slight azimuthal offset of <l> > (measured 
counter-clockwise) from the convex side of the spi- 



ral arm and an apparent gap close to the inner 
disk. The active nucleus in M51a stands out as 
a bright point source in the IR maps. Low level 
(mid-)IR emission is seen in the inter- arm region 
(about one order of magnitude fainter than the 
peaks). While the spatial coincidence between the 
CO (1-0) and hot dust emission is very good for the 
central disk, it starts to vanish for the spiral arms, 
especially when only the high surface brightness 
regions are considered. 

The polar representation (Fig. 4, middle left) 
clearly reveals gaps between the central disk and 
the spiral arms: in the northern spiral arm the 
gaps are slightly more prominent ranging from 35^' 
to 43" and in the southern arm at slightly smaller 
radii of 28^' < r < 36". After these gaps, the 
bright 24 fim spots along the spiral arms are typ- 
ically azimuthally offset from the CO emission. 
The polar cross-correlation (Fig. 5, middle left) 
shows that the offset of the hot dust emission as- 
sociated with the spiral arms starts at radii of 
r > 50^' (i.e. after the gaps) with an offset of 
~ 30° that is radially decreasing to ~ 5° and ris- 
ing again for r > 80'^ with a decline thereafter. 
A recent study by Louie et al. (2013) finds quali- 
tative similar results while the lower (13") resolu- 
tion study of Foyle et al. (2011) sees no evidence 
for offsets. This finding can be easily explained by 
the spottiness of the bright 24/im peaks that are 
clearly associated with young massive star forma- 
tion and the enhancement of 24/im emission in the 
molecular spiral arms which become mixed at the 
13" resolution that no longer resolves the spiral 
arms. It is interesting to note that also the bright 
24/im spots do not fall on a straight line suggesting 
a certain stochasticity when the initial star forma- 
tion event occurs with respect to the spiral arm 
passage. 

A low level component of the 24/im emission 
coincides with the CO emission in the spiral arms 
and extends into the downstream inter-arm region 
(Fig. 2, middle row and Fig. 4, middle left). 
The factor 2 x enhancement in the smooth compo- 
nent of the 24 /im emission in the arm compared 
to the inter-arm region is similar to the average 
arm/inter-arm difference in CO brightness. This 
suggests that the enhancement simply reflects a 
constant dust-to-gas ratio. Therefore the 24/im 
emission detected in the (molecular) arms is likely 
not tracing any deeply embedded star formation. 



Detailed modeling of the dust SED at an angu- 
lar resolution sufficient to exclude the influence of 
the bright star forming regions associated with the 
arms would be required to firmly test this inter- 
pretation. 

There appears to be only a very broad corre- 
lation between the hot dust emission traced by 
24/im and the CO emission when considering the 
entire PAWS FoV (Fig. 6, upper middle row). 
However, comparison of the three individual envi- 
ronments provides more insights. The central disk 
shows a very steep linear relation, though with sig- 
nificant scatter, while there might be an additional 
loose relation with a significantly shallower slope 
present in the spiral arms and the inter-arm re- 
gion is almost consistent with no relation at all. 
The second shallower relation is arising from the 
low level 24 /im emission seen in the 2D and polar 
representation (Figs. 2, middle left^ and 4, middle 
left). 

4.2.4- (Non-thermal) radio continuum 

The radio continuum emission at 6 and 20 cm 
in M 51 is dominated by non-thermal synchrotron 
emission (Dumas et al. 2011). Surprisingly the 
correspondence to the CO (1-0) distribution is sig- 
nificantly better than for the hot dust emission 
(Fig. 2, bottom row). While bright radio contin- 
uum emission peaks exist along the spiral arms, 
the radio continuum spiral arms coincide with the 
molecular gas spiral arms. However, the radio con- 
tinuum arms appear to be more contiguous and 
diffuse compared to the CO (1-0) emission at com- 
parable angular resolution. A wavelet decompo- 
sition of the radio continuum suggests that the 
emission arising from the spiral arms is more dif- 
fuse and not composed of many individual point 
sources (Dumas et al. 2011) confirming the visual 
impression of Fig. 2 {bottom row). 

In the polar coordinate system the brightest 
knots in the 6 and 20 cm radio continuum along 
the spiral arms (Fig. 4 middle rights only shows 
20 cm) are often offset from the gas arm and align 
with the brightest regions seen, e.g., in the hot 
dust emission traced by the 24/im map (Fig. 4 
middle left). This suggests that recent star forma- 
tion is dominating the radio excitation here while 
another process must cause the enhancement of 
the radio continuum in the spiral arms. This be- 
havior is only mildly reflected in the polar cross- 
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correlation with the CO(l-O) emission (Fig. 5, 
middle right) where the azimuthal offset is con- 
sistent with zero at most radii, except for the cen- 
ter where the radio jet is present and along the 
spiral arms for 57" < r < 80'^ Interestingly, 
the azimuthal offset is only ^ ~ 8° (compared to 
^ ~ 30° seen for the hot dust emission). The 
southern arm exhibits a bubble-like feature in the 
20 cm continuum map in the polar coordinate sys- 
tem at a radius of r ~ 65^' between an azimuthal 
angle of 220 < ^ < 240°. 

The inner disk is dominated by the well-known 
radio jet (e.g. Crane & van der Hulst 1992; Mad- 
dox et al. 2007) emanating from the central AGN, 
while the molecular ring has a radio continuum 
counterpart. Inter- arm radio continuum emission 
is detected, but at a significantly lower level. At 
6 cm a second pair of spiral arms appears mid- 
way between the prominent molecular arms. New 
EVLA data taken at 20 cm shows this feature at 
that wavelength as well (priv. comm. J. Ott). 
The location of this second pair of arms is roughly 
consistent with the location of 'magnetic arms' 
in the inter-arm region shown by Fletcher et al. 
(2011). Already Tilanus et al. (1988) noted that 
non-thermal radio emission associated with the 
spiral arms (after removal of a 'base' disk in their 
8 " resolution 6 and 20 cm short spacing corrected 
VLA data) coincides with the dust lanes rather 
than star forming giant HII complexes probed by 
Ha imaging. 

Two tentative correlations between the radio 
continuum at 6 and 20 cm and the CO emission are 
apparent in the pixel- by-pixel comparison (Fig. 6, 
two bottom rows): one steeper (present in the cen- 
ter region) and a more shallow one (in the arm 
region) which separates better in the 20cm con- 
tinuum versus CO emission plot. The situation 
in the inter-arm regions seems to be more com- 
plex, judging from the significantly larger scatter 
in the respective diagrams. The steeper correla- 
tion is also significantly stronger in the central disk 
while the shallower relation dominates in the spi- 
ral arms and (tentatively) the inter-arm region. It 
is interesting to note that in the spiral arm region 
of the 6 cm vs. CO emission (Fig. 6 middle right 
panel in second last row) there is an indication 
for a down turn at large CO fluxes which might 
point toward synchrotron losses in very dense gas. 
Better data is required to confirm this tentative 



trend. Slopes between the CO flux and 6 cm, re- 
spectively 20 cm, radio continuum flux density are 
not identical. The few points scattered along an 
almost vertical branch (esp. evident in the center 
diagram of the 6 cm emission) are likely caused by 
the prominent radio jet. 

4.2.5. PAH emission 

Our stellar continuum corrected IRAC images 
at 3.6 /im and 8/im contain two non-stellar emis- 
sion components, namely PAH emission features 
and hot dust continuum. The 3.6 /im image cov- 
ers the 3.3 /im PAH feature while the 8/im image 
captures emission from the prominent PAH fea- 
tures at 7.7 and 8.6 /im. Comparison of the stel- 
lar continuum corrected 8 /im (Fig. 3 bottom left) 
and 3.6 /im images (Fig. 3 top right) to the dis- 
tribution of the CO emission shows that the non- 
stellar 3.6 and 8/im emission is co-spatial with the 
CO emission in the ring and at smaller radii. Ne- 
glecting the nucleus (where the AGN could con- 
tribute to or even dominate the emission seen), 
very bright 3.6 (8) /im non-stellar emission is found 
in the northwestern part of the ring close to the 
brightest CO emission in the ring. As this is along 
the direction of the radio jet, it is not clear if this 
is just a coincidence. 

The spatial coincidence between non-stellar 
3.6 (8) /im and CO emission becomes less per- 
fect along the spiral arms. The brightest emis- 
sion regions at 3.6 (8) /im are just situated outside 
the molecular spiral arms. However, a roughly 
one order of magnitude fainter non-stellar 3.6 /im 
emission (~ 0.15 MJy sr~^) is associated with a 
significant portion of the molecular spiral arms. 
The fainter 8/im PAH emission in the inter-arm 
region exhibits a filamentary structure connect- 
ing the two spiral arms. We find that our CO 
emission aligns well with the brighter regions, es- 
pecially when looking at an unsharped version of 
the 8/im PAH image (Fig. 3 bottom right). The 
non-detection of CO emission from the faintest 
filaments in the inter-arm region might thus be 
due to our limited CO sensitivity. Overall, the 
spatial distribution of the PAH emission seems 
a very good predictor of the CO emission while 
deviations are seen in regions of massive star for- 
mation. 

The high surface brightness peaks seen in both 
the non-stellar 3.6 and 8/im emission along the 
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spiral arm (Fig. 4 bottom left and right) correlate 
very well with the peaks present in the 24 fim map 
(Fig. 4 middle left) indicating that in addition 
to the PAH features hot dust could contribute. 
The very similar radial profile of the polar cross- 
correlation of the non-stellar 3.6 /im emission with 
the CO emission (Fig. 5, bottom left) supports this 
notion. The cross correlation profile of the non- 
stellar 8 jam emission (Fig. 5 bottom right) exhibits 
a different shape with an azimuthal offset typically 
smaller than the one seen in the hot dust emission 
(see §4.2.3). The average offset is only ^ ^ 15° 
and shows much less pronounced variations with 
radius in the spiral arm region of r > 55'^ The 
difference in the cross-correlation profiles of the 
non-stellar 3.6 and 8 /im emission can be explained 
with different amounts of hot dust contribution. 
The distribution of the faint PAH emission is more 
obvious in the 8/im PAH emission due to better 
sensitivity (Fig. 4 bottom right). The faint PAH 
emission follows very well the gas spiral arms and 
the CO distribution in the ring. The gaps seen at 
3.6 /im in the spiral arms at roughly r ^ 35^' are 
not obvious at 8/im suggesting that the gaps are 
due to a lack of sensitivity and not necessarily a 
lack of PAH emission. It is interesting to note that 
the spiral arms are wider in the non-stellar 8/im 
emission than seen in CO, specially in the radial 
range where prominent star formation is present. 
Both non-stellar components at 3.6 and 8/im 
have two broad relations with CO emission with 
slightly differing slopes and a prominent vertical 
offset in PAH intensity (Fig. 7 two top rows). 
However, there is no obvious dividing line between 
these two relations. As shown by Meidt et al. 
(2012) the 3.6 /im non-stellar emission can be a 
combination of PAH features and hot dust, espe- 
cially, in regions of active massive star formation, 
thus a single relation is not necessarily expected. 
As the upper relation is most prominent in the 
central region (Fig. 7 middle left), this strongly 
suggests that a contribution from hot dust is at 
least adding to this relation. The ratio of 24 /im 
and non-stellar 8 /im emission (Fig. 7 lower mid- 
dle row) exhibits also a steeper slope in the center, 
consistent with enhanced hot dust contribution to 
the non-stellar 8/im light. Basically no relation 
between the ratio of 24 /im and non-stellar 8/im 
and the CO emission is seen for the spiral arms, 
suggesting that this might be indeed genuine PAH 



emission that correlates with dense molecular gas. 
As the bulk of points in the inter-arm region corre- 
sponds to higher PAH intensities than in the arm 
region, we speculate that both regimes - hot dust 
mixed with PAH as well as PAH only - are present, 
although, on a significantly lower level than ob- 
served for the center. The ratio of the non-stellar 
8 /im/3.6 /im emission shows some dependence on 
CO emission as a function of galactic environment. 
The smaller scatter of this relation in the central 
disk could be interpreted as being dominated by a 
(constant) hot dust component, while in the spiral 
arms the scatter in this ratio becomes large reflect- 
ing the large range of physical conditions sampled 
in this region. The larger scatter implies either a 
more varying hot dust contribution or an intrinsic 
variation of the ratio of the PAH features sampled 
by the 3.6 /im and 8 /im filters. 

4.3. Comparison of CO (1-0) to Star For- 
mation and Stellar Tracers 

Comparison of the morphology and distribution 
of the molecular gas which is viewed as the fuel for 
star formation (e.g. Bigiel et al. 2008, 2011; Leroy 
et al. 2008; Schruba et al. 2011) to tracers sensi- 
tive to star formation with varying age sensitivity 
is illustrative to study the relation between the 
fuel and the end-product on cloud-scales. The Ha 
emission from HIT regions traces the youngest star 
forming sites between a few and about 10 Myr (e.g. 
Whitmore et al. 2011). The thermal radio con- 
tinuum should trace the same population of stars 
responsible for the ionization of the HIT regions, 
however, it is basically unaffected by the presence 
of dust (Condon 1992; Murphy et al. 2011, 2012). 
The B and V band light from young stellar clus- 
ters probes ages out to several 10 Myr (for results 
on M51, see Chandar et al. 2011; Bastian et al. 
2005; Scheepmaker et al. 2007, 2009; Haas et al. 
2008; Hwang & Lee 2008), similarly the far- and 
near-UV emission can arise from main sequence 
stars as late as A type, again being sensitive to 
ages of up to ^lOOMyr. Finally, the H band and 
3.6 /im filters are most sensitive to a wavelength 
range where the old stellar population consisting 
of red giants dominates the spectral energy distri- 
bution (in the absence of red supergiants and/or 
AGB stars which have typical lifetimes of a few 
10 to 10^ Myrs and should also exhibit a clustered 
distribution). 
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4.3.1. HII regions 

Most of the Ha emission is closely associated 
with the molecular gas (Fig. 8, top left). In 
the center the Ha emission coincides with the CO 
emission out to a radius of r ~ 35^' (Fig. 10, top 
left). A partial gap is present adjacent to the inner 
gas disk and continues out to a radius of r ~ 41" . 
After this gap the Ha emission tends to be offset 
toward the convex side of the CO spiral arm by 
about (9 - 18)° at a given radius. It is interest- 
ing to note that the Ha emission is not forming a 
straight line but is oscillating randomly as a func- 
tion of radius within this offset range. ^ Partic- 
ularly, along the spiral arms the Ha morphology 
exhibits large diffuse envelopes and shells signal- 
ing the presence of giant HII regions (see also Lee 
et al. 2011; Gutierrez et al. 2011). These very 
bright and large HII regions tend to overlap with 
the gas spurs emanating from the gas spiral arms. 
Faint Ha emission is also detected in the inter- 
arm region with an apparent preference for being 
mid- way between the two spiral arms. 

There is no obvious relation between the Ha 
emission related to HII regions and the CO emis- 
sion when comparing it on a pixel-by-pixel basis 
in the full PAWS FoV (Fig. 12, top row). A closer 
inspection of the individual region shows that Ha 
emission in the central disk exhibits some complex 
relation, while we see basically no relation between 
Ha and CO line emission in the arm and inter- arm 
region. 

4.3.2. (Thermal) radio continuum 

Unlike the radio continuum at longer wave- 
lengths (see §4.2), the 3.6cm radio continuum is 
distributed in a more patchy fashion and domi- 
nated by bright emission peaks (Fig. 8, top right). 
This suggests that the thermal contribution in 
these regions is higher. The radio jet is still the 
brightest feature seen in the central part. The con- 
tinuum emission is co-spatial with the CO emis- 
sion in the central disk. Diffuse, faint 3.6 cm ra- 
dio emission is associated with the gas spiral arms 
(Fig. 10, top right), while the bright peaks along 
the spiral arms show a similar offset from the CO 
emission as the Ha emission. About half of these 
peaks can be directly related to the brightest HII 



^The application of the polar cross-correlation was not suc- 
cessful due to the patchiness of the Ho; emission. 



regions. The polar cross-correlation with the CO 
emission (Fig. 11, top left) gives a very similar 
radial profile to the one seen for the 20 cm con- 
tinuum with a similar azimuthal offset of ~ 10°. 
The lack of strong 3.6 cm emission within the CO 
spiral arms implies that no strong star formation 
is hidden in the molecular gas spiral arm. 

It is interesting to note that the non-thermal/thermal 
separation of their 8 " resolution 6 and 20 cm ra- 
dio imaging led Tilanus et al. (1988) to conclude 
that the thermal emission solely arises from giant 
HII complexes without any evidence for such gi- 
ant HII complexes being hidden in the dust lanes. 
They did not rule out the existence of signifi- 
cantly smaller HII regions based on the sensitivity 
achieved in their radio maps. Given the fact that 
the radio imaging used here has an rms 2.5 x and 
11 X lower while achieving a 16x and 35 x higher 
spatial resolution (in beam area) at 6 and 20 cm, 
respectively, and that no prominent bright peaks 
have appeared at our resolution inside the spiral 
arms suggests that not much hidden massive star 
formation is occurring inside the gas spiral arms. 

Similar to the 6 and 20 cm radio continuum 
emission, the 3.6 cm radio continuum shows three 
broad relations with the CO emission in the pixel- 
by-pixel diagrams (Fig. 12, 2nd row from top): 
a) a large range in radio flux density for a fixed 
low CO intensity, b) a steep relation, and c) a 
possibly very shallow or even no relation between 
(low) 3.6 cm flux density and CO intensity. The 
explanation is analogous to the 6 and 20 cm radio 
continuum: the radio jet is causing the large range 
of radio flux for a given CO intensity in the center, 
while the steeper relation that is mainly present in 
the center is related to the larger amount of star 
formation whereas the very shallow one is exclu- 
sively seen in the arm and inter- arm regions and 
can be interpreted as non-thermal emission (see 
discussion §5.2). 

4.3.3. Young stellar clusters and stars 

The presence and location of young stellar clus- 
ters relative to the molecular gas is best seen in 
the HST B and V bands (Fig. 8, middle pan- 
els). The stellar light is significantly stronger in 
the center reflecting the presence of a stellar bulge. 
Prominent dust lanes are seen in both the B and 
V band image tracing the spiral arms as well as an 
intricate web in the center. As the CO emission 
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very well follows these dust lanes, a clear anti- 
correlation with the distribution of young stellar 
clusters is seen (Fig. 10, middle left). As the cor- 
respondence of CO emission to dust lanes is also 
very good in the center, dust lanes appear to be a 
superb predictor of the location of CO emission. 

The B band emission and CO spiral arms ap- 
pear to converge and meet at a radius of r ~ 78'^ 
There is a gap in the B band emission at roughly 
32^' < r < 42" that is more pronounced along the 
southern arm while along the northern arm a few 
small stellar clusters are present. The distribution 
of stellar clusters is significantly enhanced along 
the convex sides of both gas spiral arms. The dis- 
tribution is wider and brighter along the southern 
spiral arm while more clustering of stellar clusters 
is seen along the northern spiral arm. Overall, the 
width of the arms as seen in the (young) stellar 
clusters appears significantly wider than the CO 
arms. Similarly to the Ha emission, there are a 
few prominent strings of stellar clusters mid- way 
between the spiral arms in the inter- arm regions. 

The azimuthal offset between the bright CO 
and B band emission is larger beyond a radius of 
r ^ 55" and decreasing till a radius of r ^ 80'^ 
(The negative offset for radii r < 45" is similar to 
that seen for the old stellar component and dis- 
cussed there, see §4.3.4.) The V band emission 
(Fig. 10, middle right) shows a very similar pic- 
ture to the B band and has also gaps. It is in- 
teresting to note that the V band emission is basi- 
cally anti-correlated with the 24 /am. emission along 
the northern spiral arm. The overall brightness in 
the V (and B) band falls rapidly out to a radius 
of r ~ 32^' before leveling out to the disk value. 
The polar cross-correlation analysis of the V band 
continuum with the CO emission (Fig. 11, top 
right) reveals a negative azimuthal offset for radii 
of r < 42^' that becomes a positive azimuthal off- 
set beyond r > 50^' showing a similar radial de- 
pendence as the hot dust emission (Fig. 5, middle 
left)^ however, with a much less pronounced am- 
plitude, i.e. a smaller azimuthal offset. 

Ultraviolet continuum arises from stars with 
ages less than ~100 Myrs and is thus a good proxy 
for recent star formation. The GALEX FUV and 
NUV images (Fig. 8, bottom panels) exhibit a 
bright, slightly patchy central disk with CO emis- 
sion often falling into regions of less bright UV 
emission. The extent of the central UV disk is 



larger than the CO disk. In addition two promi- 
nent spiral arms are present. Interestingly, these 
UV spiral arms do not connect to the central disk 
but show a broad gap. The UV light along the 
spiral arms is offset towards the convex side of the 
spiral arms seen in CO emission and shows sev- 
eral emission peaks that are embedded in a fainter 
smooth component that is best visible in the NUV 
continuum along the southern arm. 

As is the case for the HII regions there is no 
obvious relation between the emission from young 
stellar populations (using the HST B band as a 
proxy) and the CO emission (Fig. 12, middle row). 
The large range in B band fiuxes in the pixel-by- 
pixel diagram for the center refiects the presence of 
the stellar bulge. As the emission from the stellar 
disk and individual stellar clusters is more homo- 
geneous in intensity, the range is much lower in 
the arm and inter-arm regions. It is interesting to 
note that the anti-correlation seen in the images is 
reflected in the negative relation between HST B 
band flux density and CO flux present in the arm 
region, i.e. for a higher CO flux less B band light 
is observed. 

4-3.4- Old stellar population 

The old stellar population should be a good 
proxy for the underlying stellar potential and thus 
the possible driving force of the observed CO dis- 
tribution. The HST H band image exhibits a 
bright bulge that has roughly the extent of the 
molecular disk (Fig. 9, left). Bright stellar clus- 
ters are plentiful along the spiral arms, similar to 
their distribution in the HST B and V band im- 
ages (Fig. 9, it middle panels). Given the possibly 
large contribution from younger stellar clusters, 
the contaminant-corrected 3.6 /im map of Meidt 
et al. (2012) provides a better view of the old stel- 
lar population alone and its spatial relation to the 
molecular gas (Fig. 9, right). The nuclear bar 
(e.g Zaritsky et al. 1993) is easily identifled and 
sits inside the molecular gas ring. The molecu- 
lar gas ring, or more precisely the tightly wound 
CO spiral arms, appear to start at its tips. Both 
spiral arms are also visible in this image imply- 
ing that the gravitational potential should have at 
least one m=2 mode where the molecular gas spi- 
ral arms are. The southern spiral arm is about 10- 
20% brighter than the northern arm. Assuming a 
constant mass-to-light ratio this would imply that 
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the southern arm is about 1.1-1.2 x more massive 
than the northern arm. The peaks evident in the 
contaminant-corrected 3.6 /im image are hkely ar- 
tifacts due to imperfect correction in the regions 
of most intense star formation (Meidt et al. 2012). 

In the polar representation a shght dip in H 
band surface intensity is evident at the location of 
the CO spiral arm (Fig. 10, bottom left)^. This 
suggests that even the H band light is affected 
by extinction there. The contaminant-corrected 
3.6 /im image (Fig. 10, bottom right) exhibits 
at small radii a sinusoidal profile which reflects 
the geometry of the inner bar out to a radius of 
r ~ 17^' at an angle of, e.g., ~ 185°, corre- 
sponding to a PA (measured north through east) 
of - 147° close to the value of PA = 139° for the 
nuclear bar derived by Menendez-Delmestre et al. 
(2007). A second modulation in the radial range of 
25" < r < 50^' is prominent at, e.g., an azimuthal 
angle of 40° < < 160°. This modulation has 
a very wide peak {0 ~ 100°). The width of this 
modulation decreases significantly at radii larger 
than r ^ 50'^ The most likely explanation is that 
in addition to the spiral arms (which start beyond 
r ^ 50^') another component of the galactic po- 
tential, like a wide oval, is present in addition to 
the spiral arms. As the width of the modulation 
at the largest radii probed is significantly smaller, 
the spiral arms likely affect a small part of the disk 
at a given radius. 

As already seen for the polar cross-correlation 
of the HST B band, both the HST H band and 
non-stellar continuum corrected 3.6 /im emission 
(Fig. 11, bottom panels) exhibit a negative az- 
imuthal offset of ^ ~ 5° compared to the CO 
emission for radii of r < 42" . This negative off- 
set can be easily understood, as these radii are 
inside the corotations of the nuclear bar and the 
spiral arms (plus the potential oval) and therefore 
it is expected that the gas response is leading the 
gravitational potential of these non-axisymmetric 
components that are most prominent in the old 
stellar light. Beyond radii of r > 50^' the profile 
is consistent with no azimuthal offset. We cannot 
exclude that there is a residual effect from light 
contributed by the younger stellar population that 
leads to the almost exact zero offset along the spi- 
ral arms as the azimuthal offset decreases from the 



^The black points are artifacts caused by hot pixels. 



B band to the stellar 3.6 /im to the H band emis- 
sion. 

As expected there is, in general, no correla- 
tion between CO emission and the light of the old 
stellar population as traced by their H band and 
contaminant-corrected 3.6 /im emission (Fig. 12, 
two bottom rows). Interestingly, the central region 
shows evidence for a linear relation between CO 
flux and old stellar light, however, with a second 
plume of points that shows no correlation with 
CO flux. A possible explanation could be excess 
emission caused by young(er) stars. This seems 
unlikely as the light distribution is very smooth in 
the H band and the contaminant-corrected 3.6/im 
images (Fig. 9). A further argument against con- 
tribution from young stars is the pixel-by-pixel di- 
agram of the HST B band light vs. CO emission 
that shows no correlation at all in the center. Thus 
it seems plausible that the stellar bulge contributes 
to an enhancement of CO emission, e.g. via heat- 
ing. On the other hand, no relation between the 
old stellar population and the CO emission is de- 
tected in the two disk regions (arm, inter-arm). 

5. The molecular ISM on cloud-scales 

In the following we will place our findings in the 
context of the different galactic environments that 
are probed by the PAWS FoV. The high spatial 
resolution of 40 pc of the PAWS data allows us 
to better probe the physical origins of the (cor- 
) relations seen among the different ISM and stellar 
tracers and the CO emission. 

Broadly speaking the energy density fields that 
affect the ISM can be categorized in radiation 
fields, cosmic rays, kinetic energy fields (e.g. large- 
scale turbulence), and magnetic fields. We dis- 
cuss only the relation between the distribution of 
the molecular gas and of the cosmic rays (§5.2) as 
well as the radiation field (§5.3), i.e. from stars 
and the AGN. The role of the kinetic energy for 
organizing the molecular gas and the injection of 
(large-scale) turbulence is discussed in companion 
papers (Colombo et al. 2013; Hughes et al. 2013a; 
Meidt et al. 2013). As no high angular resolu- 
tion, high sensitivity studies of the magnetic field 
in M 51 exist, we do not investigate the importance 
of the magnetic field for determining the molecular 
gas properties. As the gravitational potential ba- 
sically organizes both the distribution of the inter- 
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stellar medium and the energy density fields that 
affect it, we start with a brief description of the dif- 
ferent galactic component of the gravitational po- 
tential present in the PAWS FoV (§5.1). This wih 
provide a guideline when we try to separate the 
different components of the energy density field (s) 
that affect the ISM. 

5.1. Gravitational potential and the ISM 
geometry 

The shape of the gravitational potential deter- 
mines the distribution of both the ISM and the 
energy density fields, and thus the different galac- 
tic environments the molecular gas resides in. The 
old stellar population can be considered the best 
tracer for the underlying stellar potential (Rix & 
Rieke 1993). Based on the distribution of the old 
stellar component (from §4.3.4) and literature re- 
sults (e.g. Rix & Rieke 1993; Zaritsky et al. 1993; 
Beckman et al. 1996; Lamers et al. 2002; Henry 
et al. 2003; Menendez-Delmestre et al. 2007) we 
identify the following components in the PAWS 
field of view (from the center outwards): 

(a) a central super-massive black hole that is ac- 
tively accreting (e.g. Maddox et al. 2007, and 
reference therein) 

(b) a bulge with a radial extent of r < 16^' 
(Fisher & Drory 2010), and 

(c) a galactic disk that contains 

(i) an inner nuclear bar with a major axis 
length of r ~ (15 — 17)'^ and an ori- 
entation of PA - 139° (Menendez- 
Delmestre et al. 2007), 

(ii) an indication of an oval with a radial 
extent of r ~ 50 '^ basically oriented 
in north-south direction, and 

(iii) a 771 = 2 spiral pattern that is con- 
sistent with a spiral density wave in 
the inner part and a second spiral pat- 
tern that is often referred to as material 
wave in the outer part (for details see 
Meidt et al. 2013). 

At all galacto-centric radii the molecular gas 
distribution follows the stellar spiral arms very 
well. Even for the section 25^' < r < 50'^ where 



we see evidence for a stellar oval the molecular 
gas is roughly resembling spiral arms. This radial 
range is also where Henry et al. (2003) and Patri- 
keev et al. (2006) found deviations from a single 
pitch angle and we see a lack of recent star forma- 
tion (see §4.2.3, 4.3.1, and 4.3.3). Further inward 
the gas distribution is reminiscent of a starburst 
ring with two tightly wound spiral arms which 
start at the tip of the inner nuclear bar. The gas 
distribution along the inner bar reflects the orien- 
tation of the nuclear bar. Taken all together this 
strongly suggests that the observed molecular gas 
distribution is heavily influenced by the galactic 
potential. In addition, kinetic energy, via non- 
circular motions or large-scale turbulence, will be 
injected in to the ISM altering its properties and 
relative composition. However, a detailed analy- 
sis of the dynamical effects is beyond the scope of 
this paper. Various aspects of how the dynamical 
environment impacts the properties of the molec- 
ular gas are studied in a number of companion 
papers on the probability distribution functions 
(PDFs) of the CO emission (Hughes et al. 2013a), 
the GMC properties themselves (Colombo et al. 
2013) and possible dynamical implications have 
been derived (Meidt et al. 2013). 

It is worth noting that the azimuthal offsets 
^ measured by the polar cross-correlation (pos- 
itive <I> for star formation tracers and negative 
^ for tracers of the underlying stellar potential) 
are consistent with gas flow that is always inside 
the corotation resonance. This is independent evi- 
dence that multiple patterns with different pattern 
speeds must be present in the PAWS FoV. 

5.2. The relativistic ISM 

Since molecular gas is thought to be the fuel for 
star formation, numerous studies have examined 
the relationship between CO emission and empir- 
ical star formation tracers in galaxies. It is well- 
established that the far-infrared emission is corre- 
lated with the CO emission on large scales within 
galaxies (e.g., Devereux & Young 1990; Young 
& Scoville 1991, and references therein), while a 
correlation between the total CO and 20 cm ra- 
dio continuum (RC) luminosities of galaxies has 
also been reported (Rickard et al. 1977; Israel & 
Rowan- Robinson 1984; Adler et al. 1991; Murgia 
et al. 2002). The correlation between the CO 
and RC emission has been shown to exist for sur- 
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face brightness as well as intensity in the Milky 
Way, and for galaxies of diverse Hubble types (e.g. 
dwarfs and ULIRGS, Adler et al. 1991; Allen 1992; 
Leroy et al. 2005; Liu & Gao 2010). Detailed stud- 
ies on a few 100 pc resolution (Murgia et al. 2005; 
Paladino et al. 2006) have shown that the rela- 
tion is linear within galaxies as well as between 
regions of individual galaxies. The robust correla- 
tion down to small spatial scales is intriguing since 
RC emission is only indirectly linked to star for- 
mation. Ionized gas surrounding young high mass 
stars produces thermal free-free radiation at radio 
frequencies, but approximately 90% of the 1.4GHz 
emission in a normal star-forming galaxy is syn- 
chrotron radiation emitted by cosmic ray electrons 
gyrating in magnetic fields. While non-thermal 
radio emission is also linked to star formation, 
since supernova remnants (SNRs) from short-lived 
high mass stars are the primary site for cosmic ray 
production in galaxies, less than 10% of the non- 
thermal radio emission is due to discrete SNRs; 
the remaining emission is from electrons that es- 
cape their parent SNRs and diffuse into the disk 
and halo, where they are accelerated by the galac- 
tic magnetic field (e.g. Lisenfeld & Volk 2000). 

In §4.2.4, we showed that the morphology of the 
RC emission closely resembles the CO distribution 
within the PAWS FoV, even down to spatial scales 
of - 100 to 200 pc (see Fig. 2). Moreover, the 
pixel-by-pixel comparison (see Fig. 6 two bottom 
rows) between the CO and RC emission in M51 
reveals two relations with slightly different slopes 
(^ 0.7 and ~ 0.5 for the correlations at 20 cm, 
and ~ 0.9 and ~ 0.7 for the correlations at 6cm), 
which appear to originate from distinct environ- 
ments within M51's inner disk. Compared to the 
arm and inter- arm regions, the radio continuum 
emission at both 20 and 6 cm in the central region 
is brighter for a given CO fiux. This relationship 
is independent of the emission due to the radio 
AGN, which is manifested in Fig. 6 as a cloud of 
points that show no relation between radio contin- 
uum and CO emission. 

Potential explanations for the RC-CO corre- 
lation fall into two broad categories: those in 
which the correlation is due to a single stellar 
population that powers both types of emission, 
and those where the correlation arises because 
the synchtrotron emissivity is closely linked to the 
dense gas distribution. Among the former class 



of models, Adler et al. (1991) argued for a di- 
rect link between CO emission and non-thermal 
radio continuum via cosmic rays produced by mas- 
sive star formation. Supernovae inject new cosmic 
ray electrons into the ISM leading to enhanced 
synchrotron emission. At the same time the ac- 
companying cosmic ray nuclei (i.e. protons) heat 
the molecular gas (for a detailed model see, e.g., 
Suchkov et al. 1993). Allen (1992) points out that 
this would imply that the observed CO brightness 
is only sensitive to the amount of excitation by 
cosmic rays and not the intrinsic gas density. As 
discussed by Murgia et al. (2005), moreover, the 
observed range of CO intensities would imply gas 
temperatures above 1000 K if this scenario is ap- 
plied. Our finding that the radio continuum in 
M51 shows no spatial offset from the CO emis- 
sion - unlike the hot dust and Ha emission, which 
presumably trace sites of massive star formation - 
also argues against this model. 

Another possibility is that the good correla- 
tion between RC and CO emission in M51's spi- 
ral arms is due to enhanced synchrotron emission 
that arises from the secondary cosmic ray electrons 
produced in the interaction of cosmic rays with the 
dense molecular gas (as proposed by, e.g, Marscher 
& Brown 1978; Murgia et al. 2005; Thompson 
et al. 2007). In this scenario, interactions be- 
tween cosmic ray protons and II2 molecules pro- 
duce (negatively) charged pions whose decay pro- 
duces secondary (cosmic ray) electrons. Thomp- 
son et al. (2007) calculate that for gas surface 
densities of Sgas ~ 0.03 — 0.3gcm~^ (i.e. ^ 
150 — 15OOM0pc~^) the timescale for cosmic ray 
protons to lose all their energy is significantly 
shorter than the diffusion timescale to leave the 
medium. As the observed H2 gas densities in M51 
are in this range, this might be a plausible ex- 
planation for the very good spatial correlation be- 
tween CO emission and radio continuum in the 
spiral arms, although the smoother appearance of 
the radio continuum would point toward some dif- 
fusion of the cosmic ray electrons from their pro- 
duction sites. While the global spectral index for 
the entire disk of NGC 5194 is 0^20-6 = 0.9 (for 
S ~ i/~^; Dumas et al. 2011), there are pro- 
nounced differences at 15^' resolution with an av- 
erage (0^20-6) of about 0.6, 0.8 and 1.0 for center, 
spiral arm and inter-arm region, based on an spec- 
tral index map presented by (Fletcher et al. 2011, 
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; their Fig. 7b). 

Since the contribution of secondary electrons to 
the total RC emission (and the ratio between the 
total RC and CO emission) depends sensitively on 
parameters such as the efficiency of energy trans- 
fer from supernovae to CRs, the ratio of primary 
electrons to nuclei in CRs, and the fraction of CO 
emission that arises from dense rather than dif- 
fuse molecular gas, it is also uncertain whether 
such a model can explain why such a tight linear 
RC-CO correlation is observed for galaxies with 
a diverse range of Hubble types. Further anal- 
ysis, e.g. estimates for the dominant CR cool- 
ing processes and timescales for the typical ISM 
conditions (gas densities, radiation field, magnetic 
field strength) in different M51 environments, and 
a detailed comparison between the RC-CO corre- 
lation and the non-thermal radio spectral index, 
would be required to test whether secondary elec- 
tron production in molecular clouds helps to ex- 
plain the close relationship between CO and RC 
emission in some regions of the PAWS field. It is 
interesting to note that attempts to detect non- 
thermal radio continuum caused by secondary CR 
electrons from Galactic clouds has not been suc- 
cessful so far (e.g. Protheroe et al. 2008, and ref- 
erences therein). 

An alternative explanation for the RC-CO cor- 
relation is that the enhanced radio continuum 
emission is due to a coupling between the gas den- 
sity p and magnetic field strength 5, via B (x p^ . 
That is, since the fiux density of the non-thermal 
radio emission 5'i.4,nth at 1.4 GHz may be written 
(e.g. Hoernes et al 1998) as 



5'l.4,nth OCUcrIcrB^ ""*\ 



(1) 



an increase in the RC emission may be linked 
to an increased magnetic field strength (5), rather 
than an increase in the number of CR electrons. 
In this equation, Iqr is the path length through 
the synchrotron-emitting region along the line-of- 
sight, and ucr is the volume density of cosmic ray 
electrons with energy between E and E + 6E, for 
an injection cosmic ray energy spectrum N{E) = 
NoE~P^ where anth = — p/2. The expression in 
Equation 1 is fairly general, as it allows for vari- 
ations in the scale-height of the synchrotron disk, 
and does not assume a relationship between B and 



As an explanation for the local RC-CO corre- 
lation, the advantage of p-/3 coupling models is 
that the physical mechanism responsible for the 
coupling between the magnetic field and gas vol- 
ume density is truly local in nature, involving a re- 
lationship between the magnetic field and charge 
density within (nearly) neutral gas. Observations 
indicate /3 ~ 0.5 ± 0.1 globally and within galac- 
tic disks, even on the scale of individual clouds 
(-100 pc, e.g. Fiebig & Guesten 1989; Berkhui- 
jsen et al. 1993; Niklas & Beck 1997). This range 
of /3 values is in good agreement with results from 
turbulent MHD simulations {/3 = 0.4 to 0.6, e.g. 
Ostriker et al. 2001; Cho & Vishniac 2000; Groves 
et al. 2003) and with equipartition between mag- 
netic and turbulent energy densities in the ISM 
{(3 = 0.5, e.g. Ko & Parker 1989). Regardless of 
y^'s exact value, B — p coupling yields 



5'l.4,nth OC UcrIcrP' 



,/3(l-anth) 



(2) 



for the correlation between the non-thermal 
1.4 GHz emission and the gas volume density. The 
gas surface density Sgas is simply the projection of 
the volume density through the disk, so for a gas 
disk with constant scale height /gas we would ex- 
pect to observe a power-law scaling with a similar 
exponent, i.e. 



^1.4, nth OC UcrIcR^sL' 
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(3) 



Since the interstellar gas in M51's inner disk is 
almost entirely molecular, we can replace Egas in 
Equation 3 with /co- The relationship between 
the surface density of CRs, Nqr = ^cr^cr^ and 
Tjgas is somewhat uncertain (see e.g. the discus- 
sion in Murgia et al. 2005). In general terms, re- 
gions with high gas surface density are linked to 
regions with high levels of star formation activ- 
ity. On one hand, high star formation rates would 
tend to increase the supply of cosmic ray electrons, 
but on the other hand may lead to higher rates of 
convective escape. Assuming that the surface den- 
sity of CRs NcR is independent of Sgas, then for 
/3 = 0.5 and Ofnth = —0.7, we obtain a slightly sub- 
linear RC-CO correlation, i.e. 5'i.4^nth ^ ^co^? i^ 
general agreement with our pixel-by-pixel correla- 
tions in Fig. 6. The different slope and vertical off- 
set between the RC-CO correlations in the center 
and disk of M51 could be due to the higher level 
of ongoing massive star formation in the central 
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region. That is, fresh cosmic ray electrons from 
recent supernovae increase Nqr (e.g. Tabatabaei 
et al. 2007), while the shallower slope is consistent 
with the shallower radio spectral index in the cen- 
tral region (which itself may be due to a higher 
rate of convective escape via a galactic wind, or 
the fact that the CR electrons are still relatively 
'warm'). The presence of a few points scattered 
along the shallower relation in the spiral arm plot 
is consistent with this interpretation, since the 
brightest radio continuum peaks coincide with the 
brightest 24 /im emission (see §4.2.4). 

The detailed analysis of the magnetic fields and 
radio continuum emission in M51 by Fletcher et al. 
(2011) found that compression of the magnetic 
field by a spiral arm shock would significantly over- 
predict the observed moderate contrast between 
arm and inter-arm region. Therefore they con- 
cluded that no strong spiral arm shock is com- 
pressing the magnetic field, consistent with the 
analysis of Meidt et al. (2013) who, based on a 
kinematic analysis of the CO emission, find no ev- 
idence for a strong dynamical shock. However, 
higher angular resolution and more sensitive ra- 
dio continuum observations as now possible with 
the up-graded Jansky Very Large Array would be 
required to further test if compression of the mag- 
netic field is indeed not important even on the 
small scales probed by our CO data. 

5.3. Radiation field 

The radiation field within the PAWS FoV of 
M51 contains - broadly speaking - three main 
sources: the interstellar radiation field (ISRF) 
dominated by the old stellar population, young ac- 
tive star forming regions with a significant energy 
increase at short wavelengths (i.e. UV continuum) 
and the AGN itself that displays a jet feature (e.g. 
Maddox et al. 2007). The spatial distribution of 
these three components varies not only radially, 
but also with respect to the geometry of the non- 
axisymmetric components (e.g. bar, spiral) of the 
gravitational potential. 

The impact of radiation onto molecular gas 
is typically described in the form of a photo- 
dominated (or -dissociation) region (PDR) that is 
nowadays defined as a region where far-UV pho- 
tons dominate the energy balance or chemistry 
of the gas (e.g. Tielens 2008). Thus PDRs are 
not only associated with HII regions, but also 



the diffuse ISM, i.e. anywhere there FUV pho- 
tons penetrate. PDRs in a broad sense are the 
transition region from the ionized to the molecu- 
lar phase of the gaseous ISM. Energetic UV pho- 
tons from nearby hot stars (e.g. within a HII re- 
gion) cause photo-ionization and -dissociation in 
the outer layer of a molecular cloud. Due to atten- 
uation of these photons by atoms, molecules and 
dust grains the following (stratified) structure ex- 
ists from outside to inside, i.e. lower to higher gas 
densities: (a) a H+/H transition zone followed by 
an atomic HI layer, and (b) the HI/H2 transition 
zone after which Hydrogen is molecular, within the 
molecular H2 gas the transition from ionized Car- 
bon to CO (CII/CI/CO) occurs at slightly higher 
gas densities. However, the structure of GMCs 
is far from being simple spheres, recent Galactic 
observations show that GMCs rather contain fil- 
aments of varying density structures (e.g Andre 
et al. 2010; Pety & Falgarone 2000) consistent with 
expectation from the ISM's turbulent nature (e.g. 
Hennebelle & Audit 2007, 2008). Therefore it is 
observationally very challenging on the size scales 
considered here to discriminate between GMCs 
that are full PDRs (like in the center of Maffei2; 
Meier & Turner 2012) or only have PDR surface 
layers (as in the center of IC 342; Meier & Turner 
2005). (A more detailed overview on the heating 
and cooling processes occurring in the molecular 
gas is provided in Appendix §A.) 

5.3.1. The old stellar population (ISRF) 

The ISRF is mainly caused by the old underly- 
ing stellar population(s). As the stellar density of 
the old population is changing from the bulge to 
the disk as evident from, e.g., the optical-to-near- 
IR radial profiles (Muhoz-Mateos et al. 2011), so 
must the incident radiation field. In the following 
we discuss the role of the ISRF component that 
is dominated by the old stellar population in the 
relationship between the molecular gas and dust 
emission (i.e. 24/im, PAH). 

Hot dust (Very small grains): As the radiation 
of the old stellar population is significantly less 
energetic than that emitted by hot OB stars it is 
not clear how much it can affect the properties 
of the CO-identified molecular clouds that typi- 
cally have extinction values of Ay > 1 (see §4.2.2 
for details). The radial variation in the old stel- 
lar population can explain the different slopes seen 
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between 24 /im emission and CO flux. The shallow 
slope seen in the spiral arms (Fig. 6 middle right 
panel in second row) can be attributed to diffuse 
24 /im 'cirrus' emission (for a detailed discussion 
see Section 5 of Leroy et al. 2012) that is caused 
by heating from the ISRF. That is we see the short 
wavelength tail of the blackbody emission of the 
cold dust component that is well-mixed with the 
molecular gas and not hot dust heated by young 
massive stars. This immediately implies that in 
the center and in certain areas of the spiral arms 
the 24 /im emission is 'over'-luminous compared 
to the CO emission. Muhoz-Mateos et al. (2009) 
studied the radial proflle of the dust emission in 
the SINGS sample, including M 51. They flnd that 
M51's surface density of the dust luminosity and 
of the derived dust mass vary independently with 
radius, with the dust in the central 2kpc being 
signiflcantly more luminous for a given mass sur- 
face density (their Fig. 7.44). Young star forming 
regions are very good candidates for causing this 
'over '-luminosity. 

However, the situation in the central disk is 
more complex. While it is very likely that the 
ongoing massive star formation (providing signif- 
icant UV photons) is the main heating source for 
the dust, and thus the enhanced 24/im emission, in 
the molecular ring, it is not obvious what is caus- 
ing the excess emission inside the molecular ring 
(best seen at 70/im in Fig. 2 middle right panel). 
This enhanced 24(70) /im emission is indicative of 
dust hotter than that seen in the spiral arms and 
indeed Mentuch Cooper et al. (2012) flnd that the 
dust temperature in the center is about 25 K, only 
slightly lower in the (dust) spiral arms and about 
3-4 K colder in the disk unaffected by recent star 
formation. The increase in dust temperature in- 
side the molecular ring could either come from an 
increase in stellar UV emission or the AGN. Un- 
fortunately, the (F)IR imaging from Herschel used 
by Mentuch Cooper et al. (2012) has not sufficient 
spatial resolution to spatially resolve the central 
region. An obvious candidate for the stellar radia- 
tion fleld would be the older stars that make up the 
bulge of M 51 that dominates over the central disk 
for radii r < 15^' (Beckman et al. 1996). Lamers 
et al. (2002) derive an age of more than 5 Gyrs 
for the prominent 'smooth, yellowish' component 
which they attribute to be the bulge. Recently 
Groves et al. (2012) showed that the old stellar 



population in the bulge of M 31 alone can explain 
the observed increase in dust temperature without 
the need for any young (< lOOMyr) stars. Simi- 
larly, Bernard et al. (2008) observe an increase in 
dust temperature at the center of the LMC where 
the density of old stellar population is enhanced 
due to the stellar bar. The relation of CO emis- 
sion and emission from the old stellar population 
(as traced by its H band and contaminant-free 
3.6/im emission; see §4.3.4) shows a remarkable 
difference for the center and the disk (arm and 
inter-arm region): The flrst region exhibits a clear 
linear relation, while the latter are basically con- 
sistent with no correlation. This could potentially 
imply that the molecular gas might experience ad- 
ditional heating by the enhanced ISRF caused by 
the higher density of old stars. 

PAH: Given the excellent correlation between 
CO and PAH emission we suggest in the following 
that the ISRF that is dominated by the old stellar 
population is the underlying cause for this correla- 
tion. Modeling of the infrared spectral energy dis- 
tribution with the Draine & Li (2007) model using 
the PACS 70 and 160 /im data in conjunction with 
the SPIRE 250 /im imaging following the method- 
ology of Aniano et al. (2012) shows that the ISRF 
in M51 is (within a factor of two) about 11 x, 7x 
and 6.5 X the Galactic value in the central, spiral 
arm and inter- arm region of the PAWS area (B. 
Groves; priv. comm.), so likely suflicient to excite 
PAH emission. 

While the exact excitation channels and condi- 
tions for PAH emission observed in the MIR spec- 
tra of nearby galaxies are not well-understood (e.g. 
Smith et al. 2007), it is well-established that sig- 
niflcant PAH emission (similarly to mid-IR contin- 
uum arising from heated small grains) is related to 
star forming regions ranging from PDRs to HII re- 
gions (see review by Tielens 2008). Photo-electric 
heating by UV photons from these star forming 
regions or UV photons present in the ISRF is 
very likely the main excitation mechanism (Tielens 
2008). Heitsch et al. (2007) investigated the pre- 
dictive power of PAH emission for cloud structure 
and its underlying density distribution using a sin- 
gle cloud illuminated by a single O star. They flnd 
that for diffuse gas where UV photons can pene- 
trate the entire cloud, a good correlation between 
PAH emission and gas density distribution is ex- 
pected. However, at higher gas densities, UV pho- 
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tons no longer illuminate the full cloud but only 
its 'rims' and thus a mis-interpretation of the un- 
derlying gas distribution will become possible and 
most likely. Based on these theoretical findings the 
correlation between PAH and CO emission will de- 
pend on both the exact density distribution of the 
(CMC) structure and the distribution of the in- 
cident radiation field. This implies that a spatial 
correlation between CO emission and PAH emis- 
sion should exist, however, the observed intensities 
might not agree, particularly in regions of intense 
star formation. 

Our observational results showing a very good 
spatial relation between CO and PAH emission 
corroborate this interpretation. Two of our find- 
ings are interesting in this context. PAH emis- 
sion is strong even in regions where no excess UV 
fiux from nearby massive star formation is present, 
like the 'gaps' in the spiral arms seen in, e.g., UV 
emission (see §4.3.3). In addition, there is tenta- 
tive evidence that the geometry of the CO distri- 
bution, especially in the inter-arm region, is very 
similar to the geometry of the PAH emission (se 
§4.2.5). However, the pixel- by-pixel comparison 
of both PAH emission features shows no predic- 
tive power in the disk, i.e. the PAH strength is 
varying a lot for a given CO fiux. Only in the cen- 
ter, a tighter correlation is present. This points 
towards a strong link between PAH emission and 
molecular clouds. Although excess UV radiation 
from slightly (up to ~ 100 Myr) older stars could 
increase the PAH emission by providing a stronger 
(local) radiation field, we find no difference be- 
tween the arm and inter-arm region when consid- 
ering the range of PAH emission probed. Thus, 
the absence of a correlation in the disk (arm, inter- 
arm) suggests that the PAH emission is not trac- 
ing the true gas density, but rather arises from 
cloud surfaces in the disk of M51 while the inci- 
dent radiation field is not very important or PAHs 
are readily destroyed. The apparent correlation in 
the center between PAH and CO emission is at 
least partially due to the contribution of hot dust 
emission to the non-stellar emission at 3.6 /im and 
8/im as the ratio of 24/im and non-stellar 8/im 
versus CO fiux exhibits a shallower slope. 

Sandstrom et al. (2010) investigated the spa- 
tial relation between PAH emission and other dust 
and gas tracers in the Small Magellanic Cloud, 
they noted that PAH emission shows the best cor- 



respondence with CO emission tracing molecular 
gas. They offered two explanations for this corre- 
lation either PAHs are formed in molecular clouds 
or the shielding offered by the dense gas environ- 
ments prevents destruction of the PAH features. 
Our good spatial correlation between PAH and 
CO emission in all environments underlines the 
close relation between PAHs and molecular gas. 
A similar conclusion has been reached for 3 clouds 
in ICIO by Wiebe et al. (2011). Taking ah this ev- 
idence together this strongly suggests that PAHs 
might indeed form in molecular gas or at least the 
dense phase of the neutral interstellar medium and 
that the ISRF is illuminating the PAHs in the 
outer shell of GMCs. Reversely, it directly im- 
plies that maps of, e.g., the 8/im PAH emission 
have a very good predictive power of the expected 
spatial distribution of the molecular line emission 
traced by CO (1-0), however, not the expected line 
strength. As a corollary this implies that 8/im 
PAH emission should have a poor predictive power 
of the recent star formation rate (at high spatial 
resolution) . It should work fine at ~kpc resolution 
where CO vs. star formation is robust. 

5.3.2. The OB stars (massive star formation re- 
gions) 

The distributions of the molecular gas and the 
sites of ongoing and/or recent star formation, as 
traced by HII regions, young stars and stellar clus- 
ters show a correlation on the scale of galactic 
structures, e.g. the spiral arms and the central 
gas ring. However, they do on average not spa- 
tially coincide when compared at ~ 1" ^ i.e. 37 pc, 
resolution (see §4.3.1 and 4.3.3). This is expected 
as the strong UV radiation from the young stars 
and their stellar winds are believed to quickly dis- 
rupt the parent cloud (e.g. Hopkins et al. 2010, 
and references wherein). In this context it is in- 
teresting to distinguish between the star forming 
sites offset from the dense molecular gas along the 
gas spiral arms and the star forming ring where re- 
gions of young stars and dense molecular gas are 
well-mixed on the 100 pc scale. In the former case 
the expectation is that one observes the proper- 
ties of a PDR not much diluted by contributions 
from the cold ISM component whereas the ring 
provides an ideal situation to obtain global prop- 
erties of PDRs mixed together with cold molecular 
gas. 
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Spiral arms: The azimuthal offset of all star 
formation tracers frora the CO spiral arms is on 
average at least ^ ~ 10 deg for radii of r > 50^' 
(i.e. past the 'gap'; see Fig. 5). This translates 
into a spatial separation of over 300 pc for radii 
beyond 1.84 kpc between the sites of massive star 
formation, and thus OB stars, and the reservoir 
of molecular gas. Although the impact radius of 
giant HII regions onto surrounding GMCs is not 
very well known, e.g. Wilson et al. (1997) quote an 
upper limit of < 120 pc based on a detailed study 
of GMCs with associated HII regions in M33, this 
immediately implies that the radiation of young 
stars is not affecting the bulk of the molecular gas 
in the spiral arms. It is interesting that emission 
from the expected dissociation products such as 
[CII] and HI is found off the dense gas arms, espe- 
cially the northern arm (see §4.2.1). Similarly to 
the [CII] (and HI) distribution, the hot dust emis- 
sion is present in distinct locations also provid- 
ing further evidence that the UV radiation of the 
young stars is impacting the ISM, as evidenced by 
the higher dust temperature, offset from the CO 
arms (see Fig. 10 of Mentuch Cooper et al. 2012). 
The spatial resolution of the infrared data is, un- 
fortunately, not high enough to derive exact source 
sizes. However, comparison to the hot dust distri- 
bution and the location of the HII regions suggests 
that the impact is fairly localized - consistent with 
the previous findings by Wilson et al. (1997). If we 
take the Ha emission as our best predictor of the 
location where radiative and mechanical energy re- 
leased by young OB stars is impacting the molecu- 
lar gas, stellar feed-back should not strongly alter 
the properties of the bulk of the molecular gas lo- 
cated in the spiral arms. However, this does not 
rule out the existence of individual regions in the 
disk where stellar feed-back has likely significantly 
impacted the molecular gas, e.g. north-east of the 
molecular ring inside the northern spiral arm. 

Star forming ring: In the molecular ring hot 
young stars and cold molecular gas are very close 
to each other, if not spatially coincident. There- 
fore it is interesting to test if this massive ongo- 
ing star formation leaves a significant imprint on 
the molecular gas. The northwestern ring seg- 
ment is brightest in CO, [CII], and 24 /im emis- 
sion suggesting that star formation might have 
the strongest impact onto the properties of the 
molecular gas. The spatial coincidence of these 



three tracers means that photo-dissociation and 
dust heating by young stars is occurring. How 
much this excess in UV radiation is affecting the 
cold molecular gas is not straight-forward to judge. 
Therefore, it is illustrative to compare the sur- 
face density of the cold molecular gas to that of 
the warm and hot H2 as derived from Spitzer/IRS 
imaging by Brunner et al. (2008). They find that 
the H2 surface density of both the warm (100- 
300 K) and hot (400-1000 K) phase is highest at 
the northwestern part of the molecular ring. Com- 
parison to the cold H2 surface density (probed by 
our CO data and derived assuming a Galactic con- 
version factor) of ~ lOOOM0pc~^ shows that the 
warm phase makes up about 1% while the hot 
phase is less than 0.02% of the cold molecular gas 
surface density. We find for the outer northern 
spiral arm in the PAWS field that is covered by 
the IRS H2 footprint (see Fig. 1 of Brunner et al. 
2008) similar values within a factor of two. This 
implies that the OB stars are not sufficient to heat 
large fractions of the molecular gas above 100 K. 
For instance, using observations from multiple CO 
transitions Weiss et al. (2001) derive an average ki- 
netic temperature of T^in > 100 K for the molec- 
ular gas associated with the central starburst in 
the nearby galaxy M82. Thus we speculate that 
most of the molecular gas in the star forming ring 
of M51 is still in the cold (< 100 K) phase despite 
the presence of strong UV heating by young stars. 
This would imply that the massive star formation 
is not significantly altering the properties of the 
cold molecular gas. However, a more detailed and 
proper analysis using multiple CO transitions is 
required to firmly test this. 

In summary, we find no strong evidence that 
recently formed massive stars are significantly im- 
pacting the bulk of the molecular gas reservoir in 
the PAWS FoV. 

5.3.3. AGN 

There is no strong evidence that the AGN is al- 
tering the properties of the cold molecular gas on 
large scales. Scoville et al. (1998) already noted 
that the nuclear molecular gas reservoir as seen 
in its CO (1-0) emission line is concentrated in a 
ring-like distribution with a size of ~ lOOpc. This 
nuclear gas concentration brightens in higher J 
transitions of the CO molecule indicating higher 
excitation temperatures (Matsushita et al. 2004). 
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Strong HCN(l-O) line emission is detected from 
this reservoir (Kohno et al. 1996; Heifer & Blitz 
1997) further indicating that the AGN is alter- 
ing the molecular gas properties here. Based on a 
detailed kinematic analysis of high angular resolu- 
tion C0(2-l) data, Matsushita et al. (2007) sug- 
gested that the molecular gas is entrained by the 
jet. 

A jet-like structure has been observed emanat- 
ing from the nucleus at X-ray and radio wave- 
lengths with an extent of about &^ to the south 
and 13^' to the north (e.g. Crane & van der Hulst 
1992; Terashima & Wilson 2001; Maddox et al. 
2007). Ford et al. (1985) suggested based on op- 
tical line ratios and their unusual width that the 
line emission is caused by shocks where the bi- 
directed jet interacts with disk material. This in- 
terpretation is consistent with the properties of 
the X-ray gas and the radio continuum that are 
both suggesting a non-stellar source for their exci- 
tation (Maddox et al. 2007). Careful comparison 
of the distribution of the CO line emission, the 
optical extinction and the 20 cm radio continuum 
(Fig. 2) shows that the CO emission is closely 
following the optical dust lanes that are reminis- 
cent of dust lanes regularly found along the lead- 
ing side of large-scale bars. However, no obvious 
spatial correlation to the jet geometry (as seen in 
its 20cm continuum) is found. We see a lack of 
CO emission at the location of the southern bub- 
ble which could be interpreted as the jet causing 
dissociation of molecular disk material. Targeted 
sensitive molecular line observations would be nec- 
essary to judge if the jet has indeed an impact on 
the molecular gas in the disk. No evidence for an 
impact of the jet is found outside a radius of ^6" 
and off the jet-features. Thus we conclude that 
the jet is basically not affecting the molecular gas 
in the central region despite its extra X-ray and 
UV emission contribution to the radiation field. 

It is worth noting that Brunner et al. (2008) 
find that the hot (400-1000 K) H2 gas is peaking at 
the nucleus. Interestingly, the nuclear distribution 
of the hot gas phase is elongated along the radio 
jet axis in the northern direction. The authors 
speculate that heating by the AGN is causing this 



excess emission. 



5.4. The relation between molecular gas 
and (cold) dust 

The excellent correspondence between the dis- 
tribution of the CO emission and the optical ex- 
tinction as traced by the I — H color implies that 
the molecular gas and cold dust are at least co- 
spatial and very likely well-mixed. This is in agree- 
ment with the large body of work on Galactic 
molecular clouds. Previous work (Bohlin et al. 
1978; Lombardi et al. 2006) found that for extinc- 
tion values below ~ (1 — 2) mag no CO emission is 
detectable (on CMC scales), in agreement with the 
expectations from models for photodissociation of 
the CO molecule (e.g. van Dishoeck & Black 1988) 
where the shielding of CO requires certain column 
densites. Studies of resolved Galactic clouds using 
extinction measures of Ay from individual stars 
with NIR colors typically find that the ^^00(1-0) 
emission has an almost linear relation to Ay for 
2 mag < Ay < (4 — 6) mag while for large extinc- 
tion values the ^^CO emission appears to saturate 
(Lombardi et al. 2006; Pineda et al. 2008; Heider- 
man et al. 2010). 

While within the PAWS FoV the morphology 
is very similar, the correlation between the CO in- 
tensity and I — H color is not constant (see Fig. 
6 top row). The slope is changing between the 
environments, most notably from center to disk. 
A change in the color of the underlying stellar 
population can be excluded to explain the differ- 
ent slopes as no radial change in color has been 
observed (e.g. I-H color by Rix & Rieke 1993) 
nor a trend in age for the (old) stellar popula- 
tion (Mentuch Cooper et al. 2012). Thus the ob- 
served change in slope can be either caused by a 
change of the properties of the obscuring mate- 
rial, i.e. the dust, or a change in the excitation 
conditions of the CO emission line. Indeed the 
dust temperature as well as the interstellar radi- 
ation field (ISRF) are higher in the center than 
in the disk (with the exception of the spiral arms, 
Mentuch Cooper et al. 2012). Pineda et al. (2008) 
also observed different slopes as a function of envi- 
ronment within one single CMC (Perseus cloud), 
while Heiderman et al. (2010) find slightly differ- 
ent dependencies between Ay and CO intensity 
for the two GMCs they analyzed. This might in- 
dicate that the relation between Ay and CO in- 
tensity depends on several parameters of the local 
environment. 
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In order to better understand and quantify 
these differences detailed modeling of the stellar 
spectral energy distribution is important. In any 
case the observed linearity between I—H color and 
CO intensity shows that the gas and dust are well- 
mixed and that there is a preferred gas-to-dust ra- 
tio within M 51 as found by Mentuch Cooper et al. 
(2012). 

Finally, we note that the correlation between 
the tracers of hot dust and the CO emission is not 
very good in every environment. Since the dust- 
to-gas ratio is basically constant around a value of 
105 over the PAWS FoV (Mentuch Cooper et al. 
2012), the simplest explanation is that the heat- 
ing source of the dust as observed at 24 /im (and 
70 /im) is varying and therefore the luminosity of 
the dust. Obvious examples are the sites of young 
star formation that can explain the azimuthal off- 
set between the molecular gas arms and the 24 /im 
(70 /im) emission along these arms. 

6. The importance of galactic environ- 
ment for molecular gas properties 

The detailed comparison of CO emission to 
other tracers of the ISM and stellar population(s) 
at cloud-scale resolution reveals that the relations 
are changing as a function of galactic environment, 
i.e. between disk and bulge, but also arm and 
inter-arm regime. It is interesting that we find no 
close or simple correlation between the presence 
of star formation and the distribution of the cold 
molecular gas. Two interesting aspects are the 
gaps observed between the central disk and the 
spiral arms, especially in the UV emission, and the 
'clumping' of stellar clusters (clearly) seen along 
the northern spiral arm with roughly equidistant 
separations. Close inspection of the locations of 
the clusters of stellar clusters show that they are 
associated with spurs in the molecular gas distri- 
bution. This will be discussed in detail in a forth- 
coming paper. 

We find that massive star formation is not al- 
ways found at or even close to the brightest CO 
flux peaks (presumably also peaks in gas density), 
i.e. the gaps in star formation tracers between 
the center and the spiral arms. As the gaps are 
also apparent in tracers of younger stellar popula- 
tions (HII regions, embedded star formation) this 
suggests that no massive star formation has been 



occurring at this particular radial range for sev- 
eral 10 Myrs. This situation is reminiscent of the 
prominent dust and gas lanes found along bars, 
that also show no strong star formation (Sheth 
et al. 2002). However, in the case of M51 strong 
shocks cannot be evoked to explain this difference. 
The paper by Meidt et al. (2013) explains this be- 
havior by pressure support of GMCs that prevents 
the collapse and thus subsequent star formation. 

Clearly the underlying gravitational potential 
is significantly affecting the molecular gas distri- 
bution. The spatial offset from the convex side 
of the spiral arms can be explained by the pres- 
ence of a spiral density wave (e.g. Vogel et al. 
1988). To some level an offset between HII re- 
gions and their birthplace, molecular gas clouds, 
is expected. Both Schruba et al. (2010) and On- 
odera et al. (2010) report that the relation be- 
tween (extinction-corrected) Ha emission and CO 
emission in M 33 breaks on scales below ^ 300 pc 
- similar to our observed angular separation be- 
tween molecular gas and HII regions. They argue 
that the time evolution of the GMCs themselves is 
adding significant scatter, as spatial offsets caused 
by a spiral density wave can be ignored in the case 
of M 33. Another point raised by Onodera et al. 
(2010) is that the Ha emission does not neces- 
sarily trace the ionizing stars but rather the shell 
of the HII regions and therefore offsets are to be 
expected. As Galactic GMCs also show a wide 
spread of more than 2 orders of magnitude in their 
star formation activity (Mooney & Solomon 1988; 
Mead et al. 1990) the non-existence of a relation 
between the Ha and CO emission at 110 pc reso- 
lution is not that surprising. However, the current 
data does not allow us to discriminate between an 
evolutionary delay or the fact that only a fraction 
of the molecular gas will form stars. 

Despite the massive star formation present in 
the PAWS FoV, we find no strong evidence for 
a general impact of this star formation onto the 
cold molecular gas as probed by the CO emis- 
sion. In particular, the relations between CO flux 
and ISM/stellar tracers are always very similar 
for the arm and inter- arm region arguing against 
an strong effect from stellar feed-back, both ra- 
diative as well as mechanical, onto the molecu- 
lar gas over the scales probed here (of about 40- 
100 pc). Alternatively, the high molecular gas sur- 
face density in M 51 could prevent a large propaga- 
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tion/penetration into the molecular gas phase. In 
order to thoroughly test such feed-back scenarios 
more detailed studies, also of other nearby galax- 
ies with high quality, high angular resolution data, 
are required. 

7. Summary and Conclusion 

We present a multi- wavelength comparison be- 
tween the molecular gas distribution and other 
tracers of the ISM and stellar components at 
cloud-scale resolution in the central 9kpc disk of 
the grand-design spiral galaxy M 51. For our study 
we utilize new ~ 1" resolution CO (1-0) imaging 
from the PAWS project in conjunction with exist- 
ing archival data from the UV to radio regime. In 
particular, we find that 

• Based on the good spatial correlation of 
low frequency radio continuum and CO line 
emission, we argue for a physical relation 
between radio continuum and CO emission 
that does not directly depend on the forma- 
tion/presence of massive stars. This imme- 
diately implies that (at least in certain galac- 
tic environments) the relation between CO 
line emission and radio continuum is more 
fundamental than between CO emission and 
IR emission. We show that the relation be- 
tween magnetic field strength and gas den- 
sity expected from equipartition can explain 
our observations. 

• The good spatial correlation between CO 
and PAH emission seen in the entire PAWS 
FoV suggests that PAH molecules are closely 
associated with molecular gas and that even 
the old stellar population is providing suf- 
ficient excitation in the inter-arm region. 
However, the poor correlation of fiux be- 
tween CO and PAH emission implies that 
PAH emission is a poor predictor of the gas 
mass and most of the PAH emission is only 
arising from the surfaces of GMCs. 

• The spatial correlation between molecular 
gas and ongoing massive star formation is 
complex and strongly depends on galactic 
environment (i.e. the underlying gravita- 
tional potential), showing coincidence, off- 
sets and even absence. For the spiral arm, 
we conclude that due to the large average 



minimum spatial separation of ~ 300 pc be- 
tween the molecular gas arms and the sites 
of massive star formation no profound im- 
pact onto the properties of the molecular gas 
is expected, consistent with other observa- 
tions presented here. The situation in the 
star forming ring is less clear. Comparison 
to Spitzer mid-IR H2 observations hints that 
the conditions are not as extreme as in the 
centers of nearby starburst galaxies. 

• The impact of the AGN onto the molecular 
gas in the disk of M51 is not strong and 
appears to be restricted to a small region 
around the nucleus. 

• The optical-near-IR color I — K provides a 
very good mapping of the distribution of the 
molecular gas and shows a good correlation 
to the CO emission, however, with a chang- 
ing slope as a function of galactic environ- 
ment. 

This study is the first in a series of papers in- 
vestigating the molecular gas at CMC resolution 
in the disk of M51 from PAWS. It clearly demon- 
strates that future studies of galactic disks with 
ALMA will provide new insights into the physics 
of the ISM as well as the star formation process. 
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A. Heating and Cooling of Molecular Gas 

Here we provide a brief summary of the heating and coohng mechanisms of molecular gas in (star forming) 
galaxies, discuss implications for the observed CO line emission and the relation between CO intensity and 
H2 mass. For the purpose of this paper we focus on the cold dense molecular gas (with n(H2)> 10^ cm~^^ 
T < 100 K) that is normally assumed to be present in GMCs. We will make use of these findings in the 
discussion of our cloud-scale relations (see §5) and environmental effects (see §6). 

A.l. Molecular Hydrogen 

There are several excitation avenues for H2 that contribute differently depending on the exact physical 
properties of the gas (e.g. density, temperature, metallicity, dust-to-gas ratio) and its surrounding (radiation 
field, cosmic ray density, shocks) which are both expected to vary across a galactic disk. 

A. 1.1. Cosmic ray heating 

The interaction between cosmic rays - mainly protons and electrons accelerated to GeV energies and 
beyond - and H2 molecules is the prime heating source for the molecular ISM (in the absence of close-by 
ongoing star formation) and proportional to the H2 density n(H2) x the cosmic ray ionization rate (cr 
(Goldsmith & Langer 1978). Cosmic ray electrons and protons (above lOMeV) photo-ionize H2 molecules 
(as these are the most abundant particles). The free electrons collide with other H2 molecules and transfer 
some of its kinetic energy. In addition, H3" is formed from the reaction of H2"+H2 -^ H3"+H. The dissociative 
recombination into H2 is a major source of heating for the case of dense molecular gas (see, e.g.. Fig. 2 of 
Maloney et al. 1996). As most cosmic ray protons and electrons can completely penetrate molecular clouds 
as evidenced by 7-ray emission from molecular clouds (e.g. Bloemen et al. 1986)^ this heating mechanism is 
expected to dominate (or even be the only one) in the most dense and cold GMCs. 

Galactic cosmic rays are thought to be accelerated by shocks from SN remnants or even Wolf-Rayet stars 
(e.g. Kotera & Olinto 2011) and have a fairly homogenous distribution throughout the MW with a galacto- 
centric decrease (e.g. Bloemen 1989). It is conceivable that the production of cosmic rays is significantly 
increased above the solar value in the presence of massive star formation. It has been shown for the nearby 
starburst galaxy NGC 253 that enhanced cosmic rays emission originating from the massive nuclear star 
formation with SFR ~ 0.1 M© yr~^ can explain the observed substantial heating of the molecular gas 
(Bradford et al. 2003; Hailey-Dunsheath et al. 2008); a similar scenario has been proposed by Papadopoulos 
(2010) to explain the CO emission in star forming ULIRGs. These regions are called cosmic ray dominated 
regions (CRDR) and might have a distinct chemical signature (Bayet et al. 2011). 

A. 1.2. II2 formation heating 

The formation of H2 on grains (e.g. Hollenbach & Salpeter 1971) provides also heating of the gas, as 
some fraction of the released binding energy goes into kinetic energy and thus subsequently heats the gas 
(Goldsmith & Langer 1978). The heating rate is proportional to the total Hydrogen density (in molecular and 
atomic form) in quadrature times the fractional abundance of atomic Hydrogen. Goldsmith & Langer (1978) 
note that if H2 destruction by cosmic rays (see above) is taken into account the maximum heating rate from 
H2 formation could be 1/3 of the cosmic ray heating rate, however, in reality the steady-state assumption 
might not be fully correct. Additionally the exact properties of dust (grain sizes, PAHs, temperature etc.) 
might affect the H2 formation rates, especially in diffuse, i.e. low density, gas (Wolfire et al. 2008). 



^The interaction between cosmic rays and (atomic and) molecular gas can be observed via the 7-ray emission that results from 
7r°-decay and Bremsstrahlung of the comic ray nuclei and electrons. 
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A. 1.3. Dust heating 

The predominant heating mechanism for dust (in molecular clouds) is via the photo-effect (Draine 1978). 
As dust couples better to the interstellar radiation field (ISFR) than the molecular gas the dust kinetic 
temperature Tdust is typically larger than the gas kinetic temperature Tgas within a molecular cloud (Maloney 
& Black 1988). Thus collisions of H2 with dust grains could increase the kinetic energy of the gas. As the 
dust heating is governed by the available radiation field, it is expected that increasing the local ISFR by 
adding an intense radiation field from, e.g. young (OB) stars or an AGN, will lead to higher dust and thus 
potentially higher gas temperature if both are well mixed and coupled. 

The energy transfer between dust and gas via collisions can be described by the energy difference of the 
H2 molecule before and after a collision with a dust grain (Burke & Hollenbach 1983). As it is possible that 
the dust temperature Tdust is higher than the gas temperature Tgas as well as the other way round, both 
heating or cooling of the gas is possible via this mechanism. The heating/cooling rate is proportional to 
n(H2)^ X AT X T^^^g with AT = Tgas - Tdust (Goldsmith 2001). 

The radiation from the diffuse UV-visible-IR ISRF will heat dust grains to temperatures higher than the 
gas temperatures observed for the dense cores of dark clouds, however, for dust grains inside the clouds this 
ISRF will be attenuated. Goldsmith (2001) calculates that in order to achieve dust temperatures of 6-10 K a 
visual extinction to the surface of the cloud of 10-12 mag is required, consistent with observations (Kramer 
et al. 1999). He notes that a further drop in dust temperature is not expected due to heating provided by 
the re-emission from dust grains themselves. Thus in molecular clouds (without star formation) the heating 
by dust is proportional to the attenuated ISFR times the H2 density n(H2) divided by the gas-to-dust ratio. 
At the same time the dust cools with T^^^^. 

Thus significant heating from dust is only expected to occur at high H2 densities n(H2) above 10^ cm~^ 
where the coupling of gas and dust via gas-dust collisions is good and can significantly reduce the gas 
temperature while only slightly increasing the dust temperature (Goldsmith 2001). Examples might be 
dense clouds close to or with embedded star forming regions (with OB stars) as such clouds have elevated 
temperatures of 40 - 80 K (e.g. Goldsmith & Langer 1978). Recent simulations confirm this picture (e.g. 
Juvela & Ysard 2011), but the authors caution that a varying grain size distribution within clouds can affect 
the gas-dust coupling and thus the resulting gas temperature by 1-2 K in dense core. 

A. 1.4- UV heating 

It is expected that the intense (UV) radiation from young stars in star forming regions heats the gas. Two 
mechanisms are generally considered: heating via the grain photo-electric effect (Tielens & de Jong 1979) 
and direct UV pumping of H2 molecules followed by collisional de-excitation (Hollenbach & McKee 1979; 
Black & Dalgarno 1976; Sternberg & Dalgarno 1989). Detailed models of these so-called Photon Dominated 
Regions (PDRs) (Tielens & Hollenbach 1985) have been derived where UV photons penetrate through the 
surface into the cloud while being attenuated and produce layers of ionized, atomic and warm molecular 
gas around the cold molecular gas. How far a UV photon can penetrate into a cloud depends on the dust 
opacity. 

Gas heating via the photo-electric effect has been shown to occur by PAH molecules and small graphite 
grains (grains with size a < lOOA; Bakes & Tielens 1994). About 50 % of the UV photon energy is absorbed 
by these small dust particles. Absorption of a far-UV photon by such a dust grain may lead to ionization and 
thus ejection of an electron. The excited electron loses energy through inelastic collisions with Carbon atoms 
inside the grain on its way to the surface and will leave the grain with some kinetic energy. Via inelastic 
collisions with gas molecules the thermal energy of the gas will be raised (e.g. Tielens & Hollenbach 1985; 
Bakes & Tielens 1994). The photo-electric effect on interstellar dust grains should dominate the heating of 
the neutral atomic ISM and play an important role in the warm inter-cloud medium and in PDR regions of 
cold molecular clouds (i.e. UV illuminated regions such as HII regions, planetary nebulae etc.). The heating 
efficiency could be up to 3% (for neutral grains; lower in grains that are charged). 
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Dust plays an important role in attenuating the UV radiation (at optical depths Ay ~ 1) that controls gas 
heating in order to prevent higher gas temperatures (Maloney & Black 1988). For a simple slab geometry of 
the molecular gas, the grain-photo-electric heating is only important in the outermost layer of a molecular 
cloud in the presence of radiation from the ISRF (dominated by old stellar populations), its contribution 
could become significant in the case of intense radiation fields from young stars (e.g. for the dense clouds in 
the center of M82, Maloney & Black 1988). At low column densities as they are present at the outer surface 
of molecular clouds photo- ionization heating can become important. A comparison of contemporary PDR 
codes summarizing the main characteristics deemed important is presented by Rollig et al. (2007). 

A. 1.5. X-ray heating 

X-ray dominated regions (XDRs) are expected to be present close to nearby AGN where the strong X-ray 
radiation from the accreting black hole is impinging surrounding molecular clouds. A recent model of XDRs 
is presented by Meijerink et al. (2006). While XDR modeling is analogous to PDR models, it is interesting 
to highlight some profound differences: the assumed input spectrum is much harder for an AGN than for 
hot, e.g. OB, stars, and the main heating source is direct photo-ionization of the gas which produces fast 
electrons. These electrons collisionally excite the H and H2 gas which emits Lyman a and Lyman- Werner 
photons that ionize atoms and ionize and dissociate molecules (e.g H2 and CO) (Maloney et al. 1996). As 
hard X-rays can penetrate deeply into dense clouds the heating efficiency in XDRs is close to unity, however, 
the X-ray heating and ionization rates depend on the actual slope of the X-ray spectrum (Meijerink et al. 
2006). 

In this context it is interesting to note that Lepp & McCray (1983) were among the first to consider 
the impact of a X-ray source inside a molecular cloud, as most stars are sources of soft X-ray emission and 
a number of accreting white dwarfs should be present in each GMC. While the contribution from X-ray 
heating via the old stellar population is hard to assess separately from cosmic ray heating, Glassgold et al. 
(2012) revisited the impact of cosmic ray and X-ray heating on cold dense clouds based on new findings of 
the H3" ion abundance in Galactic clouds. They find that the combined cosmic ray and X-ray heating (also 
referred to as chemical heating) significantly changes with the properties of the medium (molecular/electron 
fraction, total density of H nuclei) but has only a weak temperature dependence. 

A. 2. Carbon-Monoxid (CO) Line Emission 

It is generally assumed that the H2 and CO molecular gas is in thermal equilibrium so that the kinetic 
temperatures are the same. Detailed calculations for dense, fully shielded molecular regions (i.e. Ay > 10) 
assuming a steady-state chemical composition predict that the CO molecules dominate the cooling in the 
low density (n(H2) < 10^ cm~^) and low temperature (Tkm < 40 K) regime while other molecules become 
the dominant coolants at high densities and high temperatures (Neufeld et al. 1995). 

The radial distribution of H2 and CO molecules across the surface into the center of molecular clouds is 
varying as the dissociation conditions of CO and H2 are not the same, van Dishoeck & Black (1988) provide a 
detailed description and modeling of CO photo-dissociation: While the unattenuated photo-dissociation from 
stellar UV light for CO is high, substantial reductions happen inside the clouds including: a) self-shielding, 
i.e. growth toward saturation of the absorption in the line itself with increasing depth and column density, 
b) mutual shielding, i.e. the blending of lines of the isotopic species with those in the same band of the 
more abundant ^^CO molecule, c) shielding by coincident lines of H and H2, and d) attenuation of radiation 
by dust particles and Carbon atoms which significantly decreases the photo-dissociation of CO even inside 
diffuse clouds. Recently, Visser et al. (2009) up-dated the calculations from van Dishoeck & Black (1988) 
using the latest laboratory data and fully treating the heavier CO isotopologues ^^CO, C^^O, and^^C^^O. In 
addition, they extended the calculations to higher excitation temperatures as well as densities covering PDR 
and circumstellar disk environments. 

However, the exact relation between CO line emission (particular the low-J transition) and H2 gas prop- 
erties is still not very well understood. Recently, Glover et al. (2010) modeled the CO formation. They 
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found that both molecular species are expected to form quickly, however the relationship between CO abun- 
dance and gas density is more complex. The CO abundance primarily depends on the photodissociation rate 
rather than the gas density (Glover & Low 2011). Clark et al. (2012) investigate the timescales required 
for a molecular cloud to become CO detectable starting with a pure atomic phase in a turbulent medium. 
They find that the timescales depend critical on the conditions of the medium out of which the cloud forms 
with CO typically being only detectable about ~ 2 Myr before the onset of star formation. As the neutral 
gas disk of M51 covered by the PAWS observations is molecular dominated (e.g. Schuster et al. 2007) it is 
unclear how to adapt these modeling results to the specific case of M51. 

A. 3. The relation between CO line emission and Molecular Hydrogen Mass 

In order to be able to relate the observed CO line emission to the amount of H2 mass present one needs 
to understand how the molecular gas is cooling. As the type of main coolants is changing with temperature, 
gas density, and relative abundance of the coolants (which themselves depend on the ionization state of the 
gas, the intrinsic metallicity of the gas as well as other chemical processes), the observed CO emission per 
H2 is expected to vary as well (e.g. Goldsmith & Langer 1978). Typically one expects the following type of 
coolants (ignoring metallicity effects): The cold (with temperatures of 10 - 40 K or below 100 K) molecular 
gas cools via atomic and molecular line emission, however, the main cooling lines shift from atomic lines 
such as CI, 01 etc. in the low density regime (n(Il2)< lO^cm"^) to CO as the main cooling channel (for 
n(H2)< 3 X lO^cm"^) and a wide range of species including water, hydrides, molecular ions etc. for the 
highest densities where ^^CO can freeze out (Goldsmith & Langer 1978). 

In addition to the physical heating and cooling conditions of the molecular gas described above, the 
abundance of CO to H2 molecules becomes crucial when determining the relation between CO line emission 
and H2 mass present. Metallicity has the most profound effect on the abundance ratio, while other effects 
such as CO freeze-out for very high density gas or shocks could also alter the chemical composition of the 
molecular gas (Goldsmith 2001; Meijerink et al. 2010; Bayet et al. 2011). However, it is anticipated that 
metallicity has a fairly uniform effect while the others mentioned are expected to be much more localized 
across a galaxy. 

For the case of average molecular clouds, Maloney & Black (1988) showed that the conversion factor Xco 
between CO intensity and N(H2) column density is depending on the kinetic temperature T^in and velocity 
dispersion of the molecular cloud cFc which in turn is related to the H2 density such that: Xco ^ T^-^dc ~ 
T^.J^Y^n(H2) (for a slightly different dependence, see Shetty et al. 2011). This dependency can be easily 
understood as the line intensity is directly proportional to the excitation temperature (i.e. Ico ^ Tex) in 
optically thick clouds (where Tex ^ Tkin is valid, especially at the lower transitions), thus the ratio ^^^ is 
basically constant. For a gravitationally bound cloud the velocity dispersion is proportional to >/M, with 
the cloud mass M being proportional to the n(H2) density for a homogenous cloud (see Dickman et al. 1986, 
who showed that this assumption is also correct for ensembles of clouds) . 

The abundance of CO relative to H2 is another critical parameter in determining the Xco factor. If the 
metallicity is smaller the dust abundance is smaller thus the attenuation of the UV radiation with depth into 
the cloud is smaller. While H2 self-shielding will prevent a drastic reduction of the H2 zone, the destruction 
of CO will lead to an increase in cloud temperature (Maloney & Black 1988). The importance of this effect 
will depend on the column density. The implication is that a lowered CO abundance will cause the CO line 
to be no longer optically thick except at the cloud cores for clouds with similar column density but different 
CO abundance (e.g. only CO emission from the core of LMC clouds vs. the full cloud for MW clouds). For 
more details on the latest observational findings on the conversion factor in external galaxies we refer the 
reader to Sandstrom et al. (2012) and Leroy et al. (2013). 

Based on the continuous conversion factor of Xco = "^ ' — ^f -S5 derived by Narayanan et al. 

(2012) we expect a range of (1 — 3) x 10^^ cm~^ (Kkms~^)~^ for our measured Ico per pixel and assuming 
solar metallicity as observed for the PAWS FoV (Bresolin et al. 2004; Moustakas et al. 2010), i.e. Z' = 1. 
This is similar to the results of Schinnerer et al. (2010) using LVG analysis and the virial mass measurements 
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of individual GMCs by Colombo et al. (2013). 



B. Calculating angular offsets via polar cross-correlation 

In order to quantify the radial variations between two different tracers along the spiral arms, we employ a 
polar cross-correlation. This method allows for easy identification of variation in the position of the maximum 
correlation as a function of radius. 

We calculate the angular offset (j) between the molecular gas spiral arms traced in CO and ISM, stellar as 
well as star formation tracers at a second wavelength A via convolution of the two corresponding maps. All 
images have the same astrometry, angular resolution and pixel sampling (as described earlier in §3). Prior to 
the convolution, the images were deprojected from the sky plane into the galaxy plane using an inclination 
of i = 21° and a position angle of PA = 173° (as determined from analysis of the PAWS velocity field; 
Colombo et al. in prep.). Next the images were projected onto a polar coordinate system (see §4.2, e.g. 
Fig. 4) and placed into the spiral arm frame assuming a logarithmic spiral arm pitch angle of ip = 21°. 
The convolution is performed in log-polar coordinates. The spiral arm frame has been chosen for ease in 
presentation as the maximum of the auto correlation should occur very near the same azimuthal position in 
each radial bin. 

As our measure of the angular offset between the ridge of the spiral arm traced by CO and the second 
tracer at wavelength A, we adopt the intensity- weighted first moment of the cross-correlation (cc) profile 
tracing the location of maximum correlation. This is in contrast to using polynomial fits to the profile or 
to using the maximum itself. We prefer this approach for several reasons. Moment measurements made 
directly from the cross-correlation profile are less susceptible to systematic uncertainties introduced by, e.g., 
the choice of the order of or limits of polynomial fits to the profile, which typically limit the reliability of the 
measured offset (Foyle et al. 2010). The first moment, in particular, is preferred for its reduced sensitivity 
(compared to the maximum) to local variations in the fiux tracing structure at spatial scales smaller than 
the arm width (e.g. brightness peaks in CO emission tracing individual clouds) ^^. 

From the second moment we can obtain a measure of the intrinsic thickness of the arms. For arms 
with azimuthal intensity profiles that are roughly gaussian, the convolved profile has a measured width 
of cf'^^cf'co + cr^, where Geo and g\ are the intrinsic widths of the arms traced by CO emission or at 
wavelength A. We verified that the width at 95% maximum of the CO-brightness auto correlation profile (CO 
crossed with itself) corresponds well to the width estimated by eye from the morphology of CO brightness. 
We therefore measure the convolved arm width as i(;~0.64 a, where a is the 2nd moment of the profile. We 
adopt this width as the uncertainty on our measurement of (\) (see Fig. 5 and 11). 



-"^^We note that compared to the maximum, the first moment also conveys additional information regarding the symmetry of a 
given arm, or how the emission is distributed relative to the spiral ridge; asymmetric arms at one or both wavelengths results 
in a skewed profiles. This skewness can also be measured directly from the profile as the third moment. 
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Fig. 1. — The CO(l-O) line emission in the central 9kpc of M51a as observed by the PAWS project. The 
integrated intensity map is shown in a square-root scaling to emphasize the distribution of the faint emission 
{top left) and a linear scaling {top right; overlaid with contours at 40, 80, 160, and 320 Kkms~^ of the linear 
representation that is used in the subsequent figures). The location of the three region used for the pixel- 
by-pixel analysis is shown in the bottom left panel. Comparison of the CO intensity map to the THINGS HI 
robust intensity map tracing the atomic gas is presented in the bottom right panel. In the bottom left corner 
of each panel a scale bar representing 500 pc at the assumed distance of M51 and the CLEAN beam of the 
CO (1-0) data are shown. 
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Fig. 2. — Comparison of the CO(l-O) line emission in the central 9kpc of M 51a as observed by the PAWS project overlaid 
in contours (same as in Fig. 1, top right) on the Herschel [CII] line map tracing photo-dissociation (top left), the HST I — H 
color image tracing extinction {top right), the MIPS HiRes 24/im image tracing hot dust emission {middle left), the PACS 
70)Lim image tracing warm dust emission {middle right), the VLA 6cm image tracing thermal and non-thermal radio continuum 
{bottom, left), and the VLA 20cm image tracing mainly non-thermal radio continuum {bottom, right). In the bottom left corner of 
each panel a scale bar representing 500 pc and the CLEAN beam of the CO(l-O) data as well as the resolution of the respective 
dataset are shown. 
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Fig. 3. — Comparison of the molecular gas distribution of M51a to tracers of PAH emission. ^^C0(1 — 
0) intensity distribution overlaid onto the total 3.6/im emission {top left) and its non-stellar component only 
{top right), the 8/im non-stellar component mainly representing PAH emission {bottom left) and its unsharped 
masked version {bottom right) where an extra ^^C0(1 — 0) contour at 20 Kkms~^ is added. In the bottom 
left corner of each panel a scale bar representing 500 pc and the CLEAN beam of the CO(l-O) data as well 
as the resolution of the respective dataset are shown. 



39 



I, C0(1-0 



Sifiiif^^liftlll 






DO 



lOO 360 420 480 540 600 060 




80 240 300 360 420 480 540 600 



Fig. 4. — Polar representation of the ^^C0(1 — 0) line emission {top left^ and contours) and other ISM tracers: 
HST I — H color image {top right)^ MIPS 24/im {middle left)^ VLA 20cm {middle right)^ the non-stellar 
continuum at IRAC 3.6 /im and IRAC 8/im {bottom). 
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Fig. 5. — Results of the polar cross-correlation between the ^^C0(1 — 0) line emission {to^ left) and other 
ISM tracers. The panels show the radial profiles of the location of maximum correlation between the 
^^C0(1 — 0) intensity and the emission from other tracers, tracing the azimuthal offset ^ between the two: 
^^00(1 - 0) intensity {top left), HST I -H image {top right), MIPS 24 /im image {middle left), VLA 20 cm 
continuum {middle right), non-stellar continuum at IRAC 3.6 /im and 8/im {bottom). The thin lines represent 
the uncertainty defined as the width of the cross-correlation profile at 95% maximum correlation. 
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Fig. 6. — Pixel- by-pixel comparison at 3.0" resolution of the ^^C0(1 — 0) emission versus tracers of the ISM 
{from top to bottom): HST I — H color, MIPS 24/im emission, VLA 6cm radio continuum, and VLA 20cm 
radio continuum. The distribution is shown {from left to right) for the full PAWS FoV, and separately for 
the central 40", the spiral arms and the inter- arm region (as defined in Fig. 1). The density is given on a 
logarithmic color scale. 
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Fig. 7. — Pixel-by-pixel comparison at 3.0" resolution of the ^^C0(1 — 0) emission versus tracers of the 
PAH emission {from top to bottom): non-stellar emission at 3.6/im, non-stellar emission at 8/im, the ratio of 
24/im/8 Mm^^on- stellar 5 ^^^ ^^^ ratio of non-stellar emission at 8/im over 3.6/im. The distribution is shown 
{from left to right) for the full PAWS FoV, and separately for the central 40", the spiral arms and the 
inter-arm region (as defined in Fig. 1). The density is given on a logarithmic color scale. 
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Fig. 8. — Comparison of the CO(l-O) line emission of M 51a to SFR tracers. ■'^■^00(1 — 0) intensity distribution with contours 
overlaid onto the HST Ha line emission showing prominent HII regions {top left), the VLA 3.6cm image tracing thermal and 
non-thermal radio continuum (top right), the HST B band middle left) and V band (middle right) image tracing young clusters, 
the GALEX FUV (bottom, left) and NUV (bottom, right) image tracing UV emission from young (< lOOMyrs) stars. In the 
bottom left corner of each panel a scale bar representing 500 pc and the CLEAN beam of the CO(l-O) data as well as the 
resolution of the respective dataset are shown. 
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Fig. 9. — Comparison of the CO (1-0) distribution to the old stehar populations as seen in the HST NICMOS 
H band image (left) and the contamination corrected stellar light at 3.6 /im {rights from Meidt et al. 2012) 
that represents the underlying gravitational potential. In the bottom left corner of each panel a scale bar 
representing 500 pc and the CLEAN beam of the CO(l-O) data as well as the resolution of the respective 
dataset are shown. 
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Fig. 10. — Polar representation of the ^^C0(1 — 0) line emission (contours) overlaid onto SFR and stellar 
population tracers: HST Ha image(top left)^ VLA 3.6cm {top right)^ HST B and V band image {middle), 
HST H band image {bottom left), and the contaminant-corrected IRAC 3.6 /im image {bottom right). 
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Fig. 11. — Results of the polar cross-correlation between the ^^C0(1 — 0) line emission and SFR and stellar 
population tracers. The panels show the radial profiles of the location of maximum correlation between the 
^^C0(1 — 0) intensity and the emission from the other tracers, tracing the azimuthal offset ^ between the 
two: VLA 3.6cm continuum {top left), HST B band image {top right), HST H band image {bottom left), 
and the contaminant-corrected IRAC 3.6 /im image {bottom right). The thin lines represent the uncertainty 
defined as the width of the cross-correlation profile at 95% maximum correlation. 
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Fig. 12. — Pixel-by-pixel comparison at 3.0" resolution of the ^^C0(1 — 0) emission versus tracers of stellar 
light {from top to bottom): HST Ha emission, 3.6cm radio continuum, HST B band, HST H band, and 
the stellar emission at 3.6/im. The distribution is shown {from left to right) for the full PAWS FoV, and 
separately for the central 40", the spiral arms and the inter-arm region (as defined in Fig. 1). The density 
is given on a logarithmic color scale. 
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Fig. 13. — Azimuthal profiles of the polar cross-correlation between ^^C0(1 — 0) intensity and other ISM 
tracers drawn in the radial range of 30^' < r < 80'' in steps of 4.5'': ^^C0(1 - 0) intensity {top left), HST 
I — H color {top right), MIPS 24 /im image {middle left), VLA 20cm continuum {middle right), non-stellar 
continuum at IRAC 3.6 //m and 8//m {bottom). The radial position of each profile is indicated by gray-scale 
and line thickness, from small (black and thin line) to large (light gray and thick line) galacto-centric radius. 
Each profile is normalized to the maximum difference between the measured cross-correlation (cc) profile 
and the average (< cc >) signal at that radius. 
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Fig. 14. — Azimuthal profiles of the polar cross-correlation between ^^ CO (1 — 0) intensity and SFR and stellar 
population tracers drawn in the radial range of 30^' < r < 80^' in steps of 4.5'^ VLA 3.6 cm continuum 
{top left), HST B band {top right), HST H band {bottom left), and the contaminant-corrected IRAC 3.6 //m 
emission {bottom right). The radial position of each profile is indicated by gray-scale and line thickness, 
from small (black and thin line) to large (light gray and thick line) galactocentric radius. Each profile is 
normalized to the maximum difference between the measured cross-correlation (cc) profile and the average 
(< cc >) signal at that radius. 
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Table 1 
Basic Parameters for M51 (NGC5194) 



Parameter Value Comments 



Hubble type 




SA(s)bc pec 


NED 


Center 




13:29:52.7087 +47:11:42.789 


(J2000), (1) 


Vsys(LSR) 




471.7±0.3km/s 


radio convention, (2) 


Distance 




7.6 Mpc 


l"=36.8pc, (3) 


Inclination 




22±5 


(4) 


Position angle 




172±3 


(4) 


Stellar mass M^ 




3.6 X 10^° M© 


(5) 


Atomic gas mass 


Mhi 


2.8 X 10^ Mq 


(5) 


Molecular gas mass Mmoi 


6.2 X 10^ M0 


(2) 


Cold gas mass M^ 


^as 


9.0 X 10^ Mq 




fgas = Mgas/{Mg 


as + M^) 


0.2 





Note. — Basic parameters adopted for M51 (NGC5194) and used in all PAWS related publications. 
References: (1) Hagiwara (2007), (2) Shetty et al. (2007), (3) Ciardullo et al. (2002), (4) Colombo et al. (in 
prep.), (5) Leroy et al. (2008) 
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Table 2 
multi-a data of m 51 used 



A/filter/line 


Instrument 


Resolution 


Pixel scale 
Vpx 


Source 


Comments 


CO(l-O) 


PdBI+30m 


1.16x0.97 


0.3 


PAWS 


PA -73°, (1 


J band 


2MASS 


2.5 


1 


2MASS 


used as astrometric reference 


B band 


HST ACS 


0.1 


0.05 


HST Heritage 


oflFset of 0.1^' & -0.4'^ (2) 


Vband 


HST ACS 


0.1 


0.05 


HST Heritage 


offset of 0.1^' &; -0.4'^ (2) 


I band 


HST ACS 


0.1 


0.05 


HST Heritage 


offset of 0.1^' &; -0.4^^ (2) 


Ha 


HST ACS 


0.1 


0.05 


HST Heritage 


offset of 0.1^' &; -0.4^^ (2) 


H band 


HST NICMOS 


0.2 


0.05 


M. Regan 


offset 0.1^' &; -0.4^' 


i band 


SDSS 


1.4 


0.4 


DRS7 


montage of 6 tiles (see §3.7) 


Ra 


KPNO 


~1.6 


0.305 


SINGS 


SINGS, (3) 


3.6/im 


SST IRAC 


1.90^ 


0.75 


S^G 


(4) 


3.6/im 


SST IRAC 


1.90^ 


0.75 


ICA 


(5) 


4.5/im 


SST IRAC 


1.81^ 


0.75 


S^G 


(4) 


5.8/im 


SST IRAC 


2.11^ 


0.75 


SINGS 


DR5, (3) 


8.0/im 


SST IRAC 


2.82^ 


0.75 


SINGS 


DR5, (3) 


24/im 


SST MIPS 


6.43^ 


1.5 


SINGS 


DR5, offset of 0.25'' &; 0.5'', (3) 


24/im 


SST MIPS 


~2.4 


0.4 


HiRes 


offset of 0.25" & 0.5", (6) 


70/im 


SST MIPS 


18.74^ 


4.5 


SINGS 


DR5, Decl. offset of 0.375", (3) 


160/im 


SST MIPS 


38.78^ 


9.0 


SINGS 


DR5, Decl. offset of 0.375", (3) 


70/im 


HSO PACS 


5.67 


1.0 


HSO archive 


offset of -2.5" &-2.0", (7) 


160/im 


HSO PACS 


11.18^ 


1.0 


HSO archive 


offset of -2.5" &-2.0", (7) 


CII 


HSO PACS 


-11.5 


1.762 


HSO archive 


offset of 2.64" and -1.76", (8) 


FUV 


GALEX 


4.48^ 


1.5 


archive 


offset of 0.225" & 0.525", (9) 


NUV 


GALEX 


5.05^ 


1.5 


archive 


offset of 0.225" & 0.525", (9) 


3.6cm 


VLA+Eff 


2.4 


0.4 




(6) 


6cm 


VLA+Eff 


2.0 


0.4 




(6) 


20cm 


VLA+Eflf 


1.4 


0.4 




(6) 


HI 


VLA 


5.78x5.56 


1.5 


THINGS 


PA -68.0°, (10) 



^Aniano et al. (2011) 

Note. — Overview of the data used here and astrometric modification applied. See text for details on the astrometric 
corrections (§3.7) and on the source and processing of the data in Section §3. References: (1) Pety et al. (2013); (2) Mutchler 
et al. (2005); (3) Kennicutt et al. (2003); (4) Sheth et al. (2010); (5) Meidt et al. (2012); (6) Dumas et al. (2011); (7) 
Mentuch Cooper et al. (2012); (8) courtesy of K Croxah; (9) courtesy of F. Bigiel; (10) Walter et al. (2008). 
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